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Abstract 



We study a mixed heat and Schrodinger Ginzburg-Landau evolution equation on 
a bounded two-dimensional domain with an electric current applied on the boundary 
I and a pinning potential term. This is meant to model a superconductor subjected to 

an applied electric current and electromagnetic field and containing impurities. Such 
a current is expected to set the vortices in motion, while the pinning term drives 
them toward minima of the pinning potential and "pins" them there. We derive the 
limiting dynamics of a finite number of vortices in the limit of a large Ginzburg-Landau 
parameter, or e — )• 0, when the intensity of the electric current and applied magnetic 
I field on the boundary scale like |loge|. We show that the limiting velocity of the 

■ vortices is the sum of a Lorentz force, due to the current, and a pinning force. We 

. state an analogous result for a model Ginzburg-Landau equation without magnetic 

field but with forcing terms. Our proof provides a unified approach to various proofs 
of dynamics of Ginzburg-Landau vortices. 



keywords: Ginzburg-Landau, vortices, vortex dynamics, pinning, critical current 
MSG classification: 35Q99, 35B30, 35B99 



^ ' 1 Introduction 
1.1 The model 

In this paper we study the dynamics of vortices in a superconductor with applied magnetic 
field and electric current in addition to possible pinning effects, under the following mixed 
heat plus Schrodinger (or complex) flow in a bounded two-dimensional domain: 

\a + if3 \\oge\){dtU + i^u) = Aau + ^{h - \u\^) in VL 

cr((9tA + V$) = V-^/i + (ra, Vam) in 

h = Hex on dQ 

^Vau ■ ^ = iuJex ■ ^ on c?f2. 
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Here f2 is the bounded two-dimensional domain representing the region occupied by the 
superconducting sample. The unknown functions are the triple {u,A, $), where u : Q ^ C 
is the "order parameter," A : f2 — is a vector potential of the magnetic field, itself given 
by h := curl A, and $ : — t- M is the scalar potential associated to the electric field, itself 
given by E := —{dtA + V$). This is a gauge theory, i.e. {u,A,^) are only known up to 
gauge-transformations of the form u ^ ue*^, A A + V^, $ i— > $ — dt^ where ^ is smooth. 
The covariant gradient denotes V — iA. For a vector X G we write X-^ = (— X2, Xi), 
and for the perpendicular gradient we write V'^h = (Vh)-'-. Also, (•, •) denotes the scalar 
product in C defined by (a, 6) = lZ{a)lZ{h) +X(a)X(6), and {a,X) for a e C and X G 
stands for the vector in with components (a, Xi) and (a,X2). 

The function h{x) is interpreted as a pinning potential. We assume that 6 : H — t- M is a 
smooth function satisfying 

< inf h{x) < b{x) < snpb{x) < 00. (1.2) 

The situation without pinning corresponds to the case 6 = 1. 

Let us explain the meaning of the various parameters in the equation. We assume a > 0, 
cr > 0, and /3 G M. When /3 = and 6 = 1, these equations are the Gorkov-Eliashberg system 
(see |GE] ), which are the standard gauge-invariant heat flow version of the Ginzburg-Landau 
equation. The case a = 0, /3 > would correspond to a pure gauge- invariant Schrodinger 
flow. Here, for the sake of generality, we consider a > and (3 real, which corresponds to 
a mixed flow or "complex Ginzburg-Landau," also commonly considered in the modeling of 
superconductivity |Dot IKIK] . 

The parameter e > is a small parameter, equal to the inverse of k, the "Ginzburg- 
Landau parameter" in superconductivity, which is a material constant defined as the ratio 
between two characteristic length scales. We will be interested in the asymptotic limit 
£ — )■ 0, corresponding to "extreme type-II superconductors." The parameter a is called the 
conductivity. Note that the parameters a,/3,a as well as the function 6(0;) are assumed to 
be independent of e. 

The boundary conditions are what make this equation quite specific: they are meant to 
account for an applied normal current, as well as an applied magnetic field. To account for an 
incoming flow of normal current, the applied (or exterior) field H^x has to be inhomogeneous 
on the boundary. Then the incoming current /ex is given by the static Maxwell equation 

V^Hex = —lex- (1-3) 

These vector fields are defined a priori in the whole (complement of f2), but only the data 
of Hex on the boundary is needed in the equation. The data only provide the information 
of lex ■ ^, which is relevant as the normal component of current. The vector field J ex has 
the same nature as an incoming current; it can serve to model an applied voltage or surface 
charges. For a discussion of these choices of boundary conditions, we refer to Tice [Ti] , where 
they were introduced and justified. Note that the more common situation of the equation 
with applied magnetic field of intensity hex on the boundary but no "applied current" can 
be retrieved by setting Hex = hex (spatially constant) on dfl (then lex = 0) and Jex = 0. To 
simplify the dependence on e we make the structural assumptions that 

Hex = llogel H, and Jex = |loge| J (1.4) 
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for iJ, J smooth from dVt to M, respectively. In [Ti] the case where Jex = J and Hf,^ = H 
are independent of e was treated. Some larger applied fields and currents were also treated, 
but the arguments could not be extended to fields as strong as the |loge| scale that we 
consider here. As a result of the presence of these boundary conditions, the dynamics are no 
longer dissipative (or even conservative), i.e. the energy of the system can increase in time. 
We have also added the presence of the pinning weight h and the mixed flow, which were not 
considered in [Ti]. This adds more generality, but is also quite relevant for the modeling. To 
explain this let us mention more details of the physics. 

Superconductors are particular alloys that lose their resistivity when below a critical 
temperature, allowing for permanent supercurrents that circulate without loss of energy. 
However, in the presence of applied fields or currents, point vortices may appear: these can 
be seen as the zeros of the order parameter, which all carry an integer topological charge 
called degree. As e — ?■ the vortices become point-like. 

When a current is applied, it flows through the superconductor, inducing a Lorentz force 
that makes the vortices move, which in turn disrupts the flow of the permanent supercurrents. 
This is an important problem in practical applications, where a steady flow of supercurrent 
is essential. To counter this effect, a common technique is to introduce impurities in the 
material, which create "pinning sites" that pin down the vortices and prevent them from 
moving, at least when the pinning is strong enough relative to current. The impurities are 
modeled by the nonconstant function h{x), with the effect being that vortices are attracted 
to the local minima of h. One then wishes to understand at which point the current-induced 
Lorentz force is strong enough to unpin the vortices and set them in motion; such a current 
is known as the "critical current" in the physics literature. For a deeper discussion of pinning 
and critical currents, and of the physics in general, we refer to |Tnt ICEt ICHt IBFGLVj and 
the references therein. 

As we shall explain in Section II. 4[ the analysis we develop here can also be applied to 
treat the simpler model equation, which we call "Ginzburg-Landau with forcing," given by 

(a + i \\oge\ (5)dtUe = Am, + ^(1 - \ue\^) + Vh ■ Vm, + 2i \loge\ Z-Vue + feUe (1.5) 

with, say, homogeneous Neumann or Dirichlet boundary condition. Here we will assume 
that Z, /i, and are given smooth functions from Vl into M^, M, and M, respectively (note 
that despite the slightly confusing notation, Z and h have nothing to do with the currents 
or magnetic flelds of (11.11) ). In fact, the flrst step for studying (II. ip is to make a change of 
unknown functions that "removes" the boundary condition while transforming the equation 
into one similar to (II. 5p . 

We note that the well-posedness of the Cauchy problem for (II. ip as well as (11.50 can 
easily be shown by adapting the arguments of [Ti] . 

1.2 Previous work 

There have been numerous works on the dynamics of a finite number of vortices in various 
fiows for Ginzburg-Landau. Each time the goal is to derive the limiting law, as e ^ 0, for 
motion of the n vortex points, i.e. a system of n coupled ODEs, at least before the first time 
of collision of the vortices under that law. 
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The first results of tliis type were tliose of Lin |Lilj and Jerrard-Soner |JSj for the heat 
flow of Ginzburg-Landau without magnetic field (i.e. no gauge A, $): 



u 



dtu = Au + ^il-\uf), (1.6) 



e 

and with fixed Dirichlet boundary condition, then Jerrard-Colliander |JCj treated the cor- 
responding Schrodinger dynamics. This required the well-preparedness assumption Es{u) < 
Tin \loge\ + 0(1) (resp. nn \\oge\ + o(l) for Schrodinger) where Es is the Ginzburg-Landau 
energy without magnetic field and n is the initial number of vortices, all of which have degree 
±1. After accelerating time by the factor |loge|, the limiting dynamical law of vortices as 
e — )■ under the heat flow fll.61) is 

hi = -ViW{ai, . . . ,an) (1.7) 

(respectively aj- = —ViW^ai, . . . , a„) in the Schrodinger case and without accelerating time) 
where is a function of interaction between the vortices called the "renormalized energy," 
introduced in |BBHj . The case of the heat flow for the gauged equations with spatially 
homogeneous applied magnetic field was treated by Spirn in |Sp] for small fields and Sandier- 



Serfaty |SS2j for larger fields. This corresponds to setting 6 = 1, /3 = 0, H^x = h^x and 
Jex = in (11. ip above and scaling h^x with e. Complex flows started to attract attention 
recently: the dynamics were derived for the Ginzburg-Landau equation without magnetic 
field by Miot |Mij in the case of the whole plane, and Kurzke-Melcher-Moser-Spirn |KMMS] 
in a bounded domain, and for the gauged equation by Kurzke-Spirn in [KSj . 

In the case of pinning, the only complete, rigorous results are due to Lin [Li2j . who 
derived the vortex dynamics without gauge and with a different model of pinning than we 
consider. Indeed, [Li2j derives the vortex dynamics from the equation 

1 u(l — \u\'^) 

dtu = J div{bVu) + ^ J ' \ (1.8) 

For our specific pinning model there are no complete rigorous results, only partial results 
by Jian-Song |JiSo] in the case without gauge. They correctly guess the limiting dynamical 
law, and they show that if any vortices persist, then they must concentrate their energy near 
the vortex paths, which provides evidence of pinning. However, they do not prove that the 
original vortices actually persist in time or that no new vortices nucleate, and without this 
information we do not see how they can fully derive the vortex motion law. These papers 
were preceded by formal results by Chapman- Richardson ^CRJ on the motion law of a vortex 
line in three dimensions for the full magnetic model. 

There are only a few results available in the mathematics literature for the problem with 
applied current, and none consider the effect of pinning. The case lex 7^ and Jex = was 
studied in |CIH IDuH IDu2[ IDGj . where numerical and formal asymptotic results established 
evidence of current- induced Lorentz forcing in the vortex dynamics. A stability analysis of 
the normal state {u = 0) in a model with applied current and mixed Dirichlet-Neumann 
boundary conditions was performed in |Alj . For a 1-D model of a superconducting wire with 
applied current, the existence of time-periodic solutions was studied in |RSZ] . The rigorous 
study of the vortex dynamics with boundary current (either I^x and Jex) was completed for 
the first time in [Ti] . 
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Our study here also provides a relatively simple and unified approach to several of the 
situations mentioned above, combining several ingredients (pinning, applied current and 
field, mixed flow). We will present the method of proof in Section [T31 

All the studies mentioned above make some "well-preparedness" assumption, i.e. assume 
that the initial energy is not larger than nn |log£| (in some weaker or stronger form) where n 
is the initial number of vortices. The resulting dynamical laws remain valid only for as long 
as the vortices do not collide or exit the domain. A much subtler analysis is required to lift 
these assumptions and extend the dynamics past collision times. This has been done only in 
the case of the heat flow without gauge in the series of papers [BOSlt IB0S2] for dynamics 
in and |Selt ISe2] for a bounded domain. In the present paper we will not attempt such 
an analysis. 

1.3 Main result and interpretation 

Before stating the main result we need to introduce various auxiliary functions with respect 
to which we make our change of unknown functions. Let 0o be the solution to 

-crA0o + ab(j)o = in 1] 

V0O ■ = o-^^{bJ — I) ■ V on dVl. 

Note that 0o exists, is unique, and is smooth. Also, — / ■ v = —VH ■ r, for r = z/-*- the unit 
tangent on dQ. We will work in a fixed gauge in which $e = |log£:| (po- 
Define /iq : ^ — to be the solution to the PDE 

div(^)+/.o = -^V^i-V0o inn ^^^^^^ 
ho = H on dVl. 

Again, ho exists, is unique, and is smooth. Define ^ ^ — to be the solution to the PDE 

(1.11) 



A,^o = ho in Q 
,^0 = on dQ. 



Then we define Xo : ^ 
which implies that 



^0 = V^^o, (1.12) 

cm\Xo = ho, divXn = inQ , ^ 

" " (1.13) 

Xq ■ u = 0, curlXo = H on dfl. 



Finally, define ^/'q • ^ ^ to be the solution to the PDE 

A^o = div (^^^^) inn 
VV'o ■ y = J ■ y on dVt. 

The PDE is well-posed because it satisfies the necessary compatibility condition (see (12. 5p 
for a more precise statement). The functions ipo and Xq are chosen in this way so that 
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crV0o ~ V"*" curlXo = biVipo — Xq) (see Lemma [2l3] for proof). Note that all these functions 
depend only on the parameters of the problem (fe, a, (3, a, J, H) and on the domain VL, but 
not on e. The index is used to emphasize this fact. 

We transform the equations by making the change of unknown functions 

Ve = ^e-*l'°^^l'^» B, = Ae- \\oge\ Xq. (1.15) 

The energy of (i;, B) is defined by 

F,{v,B) = \ [ b\VBv\^ + ^il-\v\Y + \cnT\B\''. (1.16) 

However, the energy conditions are better phrased on the modified energy, which we will see 
only differs from F^{v, B) by a term that is o(l): 

i / , |2 I " /-I I |2n2 



n 



F,{v,B):= I -lb\VBvr + ^il-\v\y + \cnT\B\ 



+ [ \\oge\' I VV^o -Xof- VbAVb - (5 |loge|^ &0o) + / M^lyfo ■ z/. 

(1.17) 

Finally, we need to introduce the vorticity measure, or Jacobian, of the configuration. 
Following for example |SS1] . let us introduce the "space-time Jacobian" of a configuration 
(f , B) by setting 

J = d ((if, dtv)dt + {iv, dBv) + B) (1.18) 

in the language of differential forms, where ds denotes dspace + iB with B regarded as a 
one-form, dgpace is the exterior derivative with respect to the spatial variables only, and d is 
the full space-time exterior derivative. Writing J = Vi dxi A dt + V2dx2 A dt + fidxi A dx2, 
we identify the spatial Jacobian (or vorticity) 

fi{v,B) = curl {{iv,VBv) + B) (1.19) 

and the velocity vector field V = (Vi, V2) 

V{v, B) = V(zf , dtv) + dt{iv, Vbv) - dtB. (1.20) 

Since d o d = it holds that dj' = 0, which implies the continuity equation 

dtfx + curlV = 0. (1.21) 

The vorticity /i and the velocity vector field V are typically measures concentrated at the 
vortices. For example, if fj,(t) = 27r(5^(t), then V(t) = 2Tr'y-^(t)6.y(t), which shows that in reality 
V encodes the perpendicular to the actual velocity of the vortex. 
We say that the solution is well-prepared if 

fi{Ve, Bs){0) 271 J27=l d-i^aO ^OT di = ±1, 

(1.22) 

^F,(f„i?,)(o) =Er=i^K«°)+o(i). 
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where the a° are distinct points. It is known since |BBHj that each vortex carries an en- 
ergy of at least n |loge|. With pinning this easily translates into the following estimate: if 
A^) — )• 27r ^"^-j^ with di = ±1 and G distinct points, then 

n 

F,{u„ A,)>7rJ2 Ilog^l + o(|log£|). (1.23) 

i=l 

We will refer to this as the "standard F-convergence result," but we will not prove it. It can 
be easily shown by adapting results in the literature (see for example [ASS] ). Because of the 
lower bound fll.23p it is essential that we assume (11. 2p so that the minimal energetic cost of 
a vortex is 7r(inf6) |loge| > C|log£:|, which is not o(|log£:|). Without this assumption, our 
techniques would be unable to control the number of vortices. 

Notice that the well-preparedness conditions f ll.22p amount to requiring equality in f ll.23p . 
This form of well-preparedness is relatively weak in comparison to those commonly found in 
the literature, which usually require bounds by nn \loge\ + 0(1). Again, it is easy to adapt 
results in the literature to show that, given any set of n distinct points G Q, initial data 
satisfying (ll.22p can be constructed. 

Our main result on (11. ip is 

Theorem 1.1. Let {u^, As,^^) be solutions to (II. ip and choose the gauge = |loge|(/)o. 
Assume the initial data for the solutions are well-prepared in the sense that (ll.22p holds, 
where {v,B) is given by (ll.lSp . Then there exists a time > and n continuously differ- 
entiable functions ai : [0,T^,) — VL such that the initial vortices move along the trajectories 
ai (i.e. ii{ve,Bs){t) /i(t) = 27[^-di6a,(t)) and 

adi + diPaf = -2di (V-^Vo(ai) - ^oiai)) - Vlog6(ai) 
a,iO) = al ^ ■ ' 

Moreover, is the smaller of the first collision time and domain exit time of vortices under 
this law. In addition, for allt G [0,T*) we have 

e^o loge ^ 

Remark 1.2. 

1. The last relation shows that even though the energy of the system can increase, no 
significant excess energy is created. 

2. According to Lemma l27^ and equation (I1.13p . the dynamical law (ll.24p may be rewritten 
as 

'crV-^0o(ai) + ^ho{ai 



adi + d,Paf = -2di , ^ - V log 6(0^). (1.26) 

This form highlights the fact that the first term on the right may be separated into elec- 
tric and magnetic parts since (po is the electric potential (our gauge is = |log£| (po) 
and Hq is related to the external magnetic field via the boundary condition in (ll.lOp . 
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3. Our analysis may he modified in a straightforward way to handle the case of pinning 
with some smaller parameter regimes: the case of complex multiplier of the form {a + 
iPe) with /?£ = 0(|loge|), and the case of fields J^x = jexJ o,nd H^x = h^xH with 
jex = 0{\\oge\) and hex = 0(|log£:|). When one of f3e, hex, jex is smaller, i.e. of 
order o(|log£:|), the resulting dynamics are the same as those derived in Theorem \l.l\ 
by setting P = 0, H = 0, J = 0, respectively. In particular, when hex one? jex ore both 
o(|log£|) but = P |log£:|, we get 

ahi + diPaf = -Vlog6(ai). (1.27) 

In the dynamical law we identify — Vlog6(aj) as the pinning force, which is identical 
to that conjectured in |CRl IJiSoj and essentially the same as that found with a different 
pinning model in |Li2j . The other term, —2di (V^ipQ^ai) — ^^^(aj)) := —2diZ^{ai), may be 
identified as the current-induced Lorentz force. Setting 6 = 1, /3 = 0, and employing the 
definitions (11.91) and fll.lOp in the form of the dynamical law fll.26p shows that Z is the same 
as the current-induced force identified in [Ti]. If in addition Hex = hex and J = 0, then Z is 
the same as in [SS2]. It is a simple matter to rewrite ( 11.24p as 

a. = {-2diZ^ia,) - Vlog6(a,)) - {2Z{a,) - rf,V^ log 6(a.)) , (1.28) 

which reveals an interesting feature of the dynamics: one part of the forcing depends on the 
degree (pushing d = ±1 vortices in opposite directions), and one part does not (pushing all 
vortices in the same direction). 

From fll.27p it is clear that at least one of jex and hex must be of the order 0(|log£:|) in 
order for anything other than pinning to drive the dynamics. In this way we identify the 
order of the critical current as 0(|log£:|). A more precise identification of the critical current 
is somewhat more delicate since it relies on the exact form of the applied currents, J and /. 
One possible way to make this identification is to further assume that J = XJq and / = A/q 
for some fixed Jq, Jq with A G (0, oo). The dynamical law fll.24p is then rewritten as 

adi + diPdj- = -2diX {V^M^i) - X^{a,)) - Vlog6(ai), (1.29) 

where ipo and Xq are determined by replacing J, I with Jq, Jq in the PDEs (I1.9I) - (I1.14I) . We 
can then define the critical current to be Aq |loge|, where 

Aq := inf{A > I no solution to fll.29p remains confined near the local minima of 6 }. 

(1.30) 

Clearly, Aq depends on Jo, /o, and b, and is difficult to compute explicitly because of (11. 9p - 
(I1.14P unless Jq, Iq, b, and Q have some simple forms. 

In the regime we have chosen, where the applied field and currents have strength |log all 
the forcing terms in the dynamical law (11.241) have equivalent strength, and the interaction 
between the vortices (the renormalized energy found in all the works above) is negligible 
compared to them. Hence, in contrast to the dynamics derived for weak (i.e. 0(1)) fields 
in [Ti], the inter- vortex interaction disappears for strong fields. Also, these forces are strong 
enough to make the vortices move at finite speed, without the need to accelerate time. In a 
way this makes the analysis simpler as we shall discuss below. 
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1.4 Method of the proof and case of the Ginzburg-Landau equa- 
tion without magnetic field 

As in [Ti], the first step is to transform tlie equation via tlie cliange of unknown functions 
f ll.lSp . The function ipo and the vector- field Xq are in fact chosen in such a way that the 
apphed field and current disappear from the boundary condition, but appear as additional 
terms in the right-hand of the PDEs. Note also that here the pinning has the effect of 

I 1 2 

penalizing b — \us\ , so with this change we expect to be close to 1. For the sake of 
simplicity, and to extract the relevant ideas, we now present a sketch of our proof in the 
case = 0, $e = 0, i.e. for the equation without gauge, which can also have an interest 
in itself. The change of unknown function essentially transforms the equation fll.ip for 
into an equation for Ve of the form (11. 5p (for a more accurate equation, see Lemma [27i|) with 
homogeneous Neumann boundary condition, for some vector field Z{x), some function h{x) 
(in fact h = log 6), independent of e, and some function feix) which may depend on e but 
blows up much slower than, say, l/e (this is not the optimal condition). The question is then 
to understand the effect of the forcing terms V/i ■ Vm^ and 2i \loge\ Z ■ Vu^ in this simple 
"forced" Ginzburg-Landau equation and to show the f^Us term has no influence. 

In this simpler setting the vorticity fi{u) and velocity V{u) can be recomputed from 
i^M-<^M as 

fi{u) = curl(zM, Vm) = 2{idiu,d2u) (1-31) 

and 

V{u) = V(m, dtu) - dt{iu, Vu) = 2{dtu, iVu). (1.32) 

For simplicity we will denote /x^ for yu(ue) and for V{ue). 

The classical starting point to obtain the limiting dynamical law is to look at the evolution 
of the local energy density 

e,{u) = -\Vu\' + ^ , (1.33) 

use the stress-energy tensor, and try to take the limit. However, here it is much more 
convenient to use a weighted energy density 

Since e'^ appears frequently, we set 6 = so that h = log 6, which is consistent with our 
notation for pinning. 

An easy computation gives that for any function Ue, 

dtCeiUe) = div{dtUe, Vm^) - (dtUe, AUe + ^(1 - l^el^)) • (1-35) 

If Ue is a solution to (11. 5p . we then have 

+ \\oge\ {dtUe, 2Z ■ tVu,) + feiUe, OtUe). (1.36) 
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The last term in the right-hand side can be recognized as 

-fe] (1.37) 



d ( \u^f — 1 



dt 

and the one before last as |loge| Vg ■ Z. We then deduce that 

dtUue) = div {b{dtue, Vue)) - ab\dtUs\^ + |loge| bVs ■ Z. (1.38) 
We next define the weighted stress-energy tensor associated to by 



1 — \Uf 



|2 



T, = b\Vu,(S) Wu, - e,{u,) + /2X2 , (1.39) 



and a computation yields 
divT, = 6(An, + ^(1 - \ue\^) + Vh ■ Vu, + Vu,) - ee{ue)Vh + &V/, '"'' ~ \ (1.40) 



t^F.. , _ „ _ X ^ , Mf -1 



Note that is a symmetric 2x2 tensor, so divT^ denotes the vector whose coordinates are 
the divergence of the rows of T^. If Ug is a solution to (11. 5p we deduce 



divTe = ba{dtUg, VmJ + b/3 \\oge\ {idtu^, Vm^) - |loge| b{2iZ ■ Vm^, Vm^) - ee{ug)Vh 

+ bVfe 



\Ur\ - 1 



1 I?/ P — 1 

= 6a(9iM„ Vm,) - -/3 |log£| - - e,(M,)V/i + 6V/,^^ (1.41) 

where we used f ll.3ip . The limiting law is then derived from the two relations f ll.SSp and 
f ll.4ip . which can be combined to get 

dtisiue) = - div div T, + ^ |log£| div(6K) + div(6ZV.) + - dw{e,{us)Vh) 

a la. a a 

(h U, |2 _ 1 \ 
-V/,^^ J - ab\dtUs\^ + \\oge\ bVe ■ Z. (1.42) 

To be able to take limits in these relations, a priori bounds on all terms, in particular 
energy bounds are needed. In typical situations (heat or mixed fiows) the total energy 
/ ee{us) or J ee{ue) decreases in time, but this is not the case here because of the presence 
of the forcing terms which can bring in energy (this is observable in fll.38p ). However the 
"product estimate" of |SSlj provides control of V: it tells us that 



t2 



V -x 



1 



£-s>0 I log £ I 



<liniinf-^( / I IdtUef] (1.43) 



u Jn 



for any smooth vector field X, where V is the limit of as e ^ 0. This estimate, inserted in 
f ll.38p , gives control of the growth of the energy as in [Ti] . While in the timescale of the 
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estimate shows the energy increases by at most a constant in (small) finite time, here we only 
get that it increases by at most rj |log£| in the (small) time interval [0,T^,]. This is due to the 
strength of the forcing terms (now in |log£| instead of constant) and pinning terms, and it is a 
significant difference: while the rj can be taken small enough (< vrinf b) so that no additional 
vortex can appear, there can still be an energy of the same order as the vortex energy floating 
around. The usual proofs of dynamics mentioned above |Lil[ ILi2[ iJSt Sp lTi| IMit IKMMS] 



use the fact that from a priori bounds, the energy density [^^^^^(mj) can only concentrate 
at the vortex locations and converges in the sense of measures to vr^(5ai(t) where the ai(t) 
are the vortex centers. Here we cannot use this and have to accept the possibility of another 
term in the limiting energy measure, which is not necessarily concentrated at points, and 
with which we work until we eventually can prove it is zero by a Gronwall argument. So 
we denote i/(t) the hmit in the sense of measures of the energy density ^^^^^(^^(t)); all we 
know (from the F convergence result f ll.23p ) is that if there are n limiting vortices located 
at ai{t) then 

n 

z/(t) >7r5^6(a,(t)). (1.44) 

i=l 

Returning now to the discussion of taking limits in fll.38p and fll.4ip . we note that in 
situations where there is bounded order forcing or no forcing, time must be accelerated for 
vortex motion, and one takes the limit of dt^^^ on the one hand, and equates it with that 
of divTe on the other hand. The stress-energy tensor itself is not bounded as e — >■ 0, so 
it does not have a true limit. However, divT^ has a limit in "flnite parts" (see |SS3] Chap. 
13 for this viewpoint) and one can show it is equal to the gradient of the "renormalized 
energy" of |BBHj . This is the main force driving the dynamics (possibly supplemented with 
the forcing of Z) in the dynamics derived in (TI]. It is in particular this computation in flnite 
parts that makes the proofs [ml [Ll2l |JSl |Sp[ [HI [Mil IKMMS] dehcate. 



Our present situation is somewhat easier in that sense. Because of the strong forcing, 
no rescaling in time is necessary, and in f ll.38p and (ll.4ip it suffices to pass to the limit 
in ji^\(^e{ue) and jj^^divTe simultaneously. Now -p^^ is easily seen to be bounded since 
Tg is bounded by the energy density, so it immediately has a weak limit, T. The terms in 
l^el^ — 1 go to zero in the limit by the energy control and the assumption ||/e||ci ^ ~- The 
other terms in f ll.38p and fll.4ip all have limits when normalized by |log£|, thanks to the 
energy bound and the compactness results on /i^ and provided, for example, by |SSlj . 
In particular -> /i = 27i Y^"^^^ di6a,{t) while ^ V = '27rYl'i=i^i^i^Mi) where the ai{t) 
are the (continuous) vortex trajectories, and where the number of vortices and their degrees 
remain constant on [0,T^,]. We may then assume, up to extraction, that 

Ts divT, b{dtUe,Vus) abldtUsl"^ 

— )• T, j- — )• divT, ■ ¥ p, and — ■ > ( (1-45) 



|loge| |log£:| |loge| |loge| 

in the weak sense of measures in x [0,T^,]. Dividing by |log£| in ( I1.38P and (ll.4ip and 
taking the limit as e — )■ 0, we thus find 

dtu = divp- C + bV ■ Z (1.46) 

and 

div T = ap-^(5bV -bfiZ^ -uVh. (1.47) 
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The next ingredient is to perform a Lebesgue decomposition of the various measures with 
respect to the vorticity measure /i = 27r J2^=i di5ai(t)dt. This yields 

n n 

V = uo{t)dt + ^ Ui{t)Sa,{t)dt, T = To + ^ Ti{t)Sa,{t)dt, 

i=l i=l 
n n 

divT = So + ^^Si{t)5a,{t)dt, p = po + '^Pi{t)5a,{t)dt, and (1.48) 

i=l i=l 

n 

C = C0 + 5Z0W<^a,WCit, 



i=l 



where i/i{t) ,Ti{t) , Si(t),pi(t), Q{t) are functions of time and i'Q{t)dt, Tq, So,po, Co are singular 
with respect to fi. These quantities are interrelated by fll.46p and f ll.47p . and indeed from 
the decomposition fll.47p we have that 

^0 = apo — u^Vh and Si = api — f3'Kh{ai)di(xf — 27ib{ai)diZ-^{ai) — ViVh{ai). (1.49) 

The latter equation in principle gives the velocity of the vortices, hence the dynamical law; 
however, in our situation we only know from f ll.44p that z/j > 7ib{ai), but we do not know 
the precise values of z/j, nor do we know Pi or S^. 

To determine Pi we appeal to the equation fll.46p . If we formally plug the decompositions 
fll.48p into f ll.46p and throw away all but the concentrated parts, we find that 

(71 \ / 71 \ n 71 

^z/AJ =div K^MaJ -5^CA,:+27r^6(a,)rf,a^^(a,)5,,. (1.50) 
i=l / \i=l J i=l i=l 

We then observe that 

(n \ / " \ " 

^J^iSa, I = -div I ^UiCtiSa, I + ^dtiyi6a, (1.51) 

i=l / \i=l / i=l 

in the sense of distributions. Using this in fll.50p and equating the terms that are divergences 
of Dirac deltas, we find that 

Pi = —Vihi for i = 1, . . . , (1-52) 

This formal calculation can be made rigorous by integrating fll.46p against an appropriately 
chosen test function with support that moves with the vortices, and indeed we can prove 
that fll.52p actually holds. 

We can also determine that Tj = for each z = 1, . . . , ra. If we use the decompositions 
ffTlg]) in ffTirp . then we see that 

n n 

div To + div {T^5a,) = So + Y^ SAr (1.53) 

i=l 1=1 

Testing this against appropriately chosen vector fields then allows us to deduce that Tj = 
for i = 1, . . . ,n. 
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The next step is to prove that Ui = 7r6(ai), z/q = 0, and Co = by a Gronwall argument, 
the key to which is a pair of estimates for Q and Si. The first estimate, 



t " 



1.54) 



is a simple corollary of the product estimate f ll.43p . The second estimate 

n 

y^^Sj- ai< 



Jo In Jo In 



;i.55) 



may be derived by relating the pairings divT ■ S and T : DS (here A : B = Ylik^ik^ik) 
through the divergence theorem and using an appropriate test vector field S. We integrate 
the relation fll.46p over Q and between time and t to find 



Kt) - / m + / c = / 

n Jn Jo Jn Jo 



^2'Kdib{ai)Z{ai) ■ d- 



;i.56) 



i=l 



Into this equation we insert fll.49p . (11.520 . (11.540 . the well-preparedness assumption (implying 

It 



1/(0) = 7rX;r=i ^(ai(0))), and the relation vr X;r=i = ^Ya=i &(ai(0)) + Jq T^YA=i'^^'^i'^ 



after rearranging a little we are led to 



vr262(a,) 



Jo 7o 
Using the fact that Vi > Tih{ai) in conjunction with (11.55^ . we find that 

/ Mt) + - vr6(a,(t))) + \ j j Co 

Jn -^^ ^ Jo Jn 



< 



^ c M + ^ id,n ^^(''^ - ^^(^*)) + 



:i.58l 



A simple application of Gronwall's lemma, using the fact that G i?^, then allows us to 
conclude that Co = and Uiit) = nb{ai(t)), z/o(t) = for all time, i.e. we find a posteriori 
that no significant excess energy develops and that C only concentrates. 

In the final step, we use the fact that 1^0 = and Co = to show that Tq = and 50 = 0. 
We can then return to (I1.48P to compare 

/ n \ n n 

= div To + Y,Ut)5,^^t)dt = divT = 5o + Y,S^{t)5a,it)dt = Y,S^{t)5,,(,t)dt. (1.59) 



1=1 



1=1 
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Hence Si = a.e. for all i = 1, . . . , n. Since we now know z/j, pi, and Si, we may return to 
fll.49p to deduce the dynamical law 



adi + di/3dj' = —2diZ^(ai) — Vh{ai). (1.60) 

We find that in the original equation (11.51) the forcing term V/i ■ Vm^ translates into a force 
— V/i(ai) = — Vlog6(aj) acting on each vortex and pushing it toward the local minima of b, 
while the forcing term 2i |log£| Z ■ Vug translates into a Lorentz-type force —2diZ-^{ai). We 
can write this result as 

Theorem 1.3. Let Ug be a solution to (11.51) where Z{x) is a smooth vector field, h{x) a 
smooth function and fe^x) a function satisfying WfeWc^ < ^ and assume the well-preparedness 
condition 

fi{ue){0) ^ 27r Y,"^^^ di6a^(Q) for di = ±l, 

In ee{ue){0) = vr ^"^^ e'^iai) \\oge\ + o(|loge|), 

where the aj(0) G Q are distinct points. Then there exist n continuously differentiable func- 
tions tti : [0,T^,) — 7- Vt such that the vortices move along the trajectories ai (i.e. fi{t) = 
2t^ Ylii^i^a,(t)) solving 

adi + diPdf = -2d^Z^{a{) - Vh{ai). (1.62) 

Moreover, T* is the smallest of the first collision time and the first exit time of vortices under 
this law. 

Note that a simple adaptation of this allows for the treatment of the case of the Ginzburg- 
Landau equation with pinning 

{a + i/3 \loge\)dtUe = Au^ + ^(6 - \ugf), (1.63) 



where min6 > 0, since the change of unknown function Vg = ^ transforms it into 

b A \f 

^(1 - \v,\^) + Vlog6 ■ ^v, + 



a + 2/3 \\oge\)dtVe = Av^ + ^(1 - \vef) + V\ogb -Vve + Vs^^ (1.64) 



and yields the dynamical law 

adi + di(3df = — Vlog6(aj), (1.65) 

which is identical to (ll.27p . 

In more usual settings of order 1 forcing terms and logarithmic time rescaling, the method 
of proof above can also be applied, modulo the computation of the limit of div in finite 
parts, and it thus gives a unified approach to derive the dynamical law for the heat flow, 
mixed heat plus Schrodinger flow, with or without pinning, with or without an "applied 
current" type forcing term. It avoids having to choose "clever" test-functions, as done in 
particular for mixed flows, and allows for the possibility of excess-energy developping in time. 

We note that in our analysis it is crucial that the currents or forcing terms are of strength 
I log el , but again this scaling is particularly relevant since it is precisely that for which the 
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pinning force and the electromagnetic forces are of the same order. The case of stronger 
currents raises new difficulties and is still an open question. 

The paper is organized as follows. In Section [2] we present the change of functions and 
choice of auxiliary functions which allow to transform the equation (11. ip into one resembling 
(11. 5p . Section |3] contains various estimates controlling 1 — If^l and that ensure their 
compactness; it can be skipped in a first reading. In Section H] we employ the method of [Ti] 
to show that, while the energy does not decrease, it cannot increase too quickly. Sections O 
and E] contain the core of the proof of Theorem ll.lt it is there that we derive the convergence 
results and dynamical law, roughly following the method outlined in Section 11.41 



2 Reformulating the equations 

In this section we give details on the change of functions that serve to transform the equations. 
The idea follows [Ti], but the pinning term complicates the choice of functions. 



2.1 Equations 

Studying the triple (m^, A^, ^e) is not convenient because of the pinning and boundary terms. 
We reformulate the equations to remove the appearance of the applied fields from the bound- 
ary conditions. 

Lemma 2.1. Suppose ip : Q ^ M. and X : — )■ are both smooth and satisfy the boundary 
conditions 

curlX = H on d^l 

X-u = ondQ (2.1) 

Vip ■ u = J ■ u on dVt. 

Letve = Mee-^l'°s^l'^/v^ and B, = A,-\\oge\X. Wnte Z = Vip-X and Z, = \\oge\ {V^- 
X). Then {ve, B^,^ir) solve 

{a + i \\oge\ ^){dtVe + t^eVe) = ^B^V, + ^(1 - \Ve\^) + V log 6 ■ V b^v^ + 2iZ, ■ Vb^v^ 
a{dtB, + V^e) = curl B, + b{iVe,VB,Ve) + |loge| curlX + (2.3) 



in Q, along with the boundary conditions 



curl Bs = on dVl 

V_B^f£ ■ = — ^fV logfo ■ on dVl. 



(2.4) 



Proof. The PDEs (12.21) and (12. 3 p follow directly from the equations (II. ip and the definitions 
of Bgr. The boundary conditions (12. 4 p follow from (12. ip and the boundary conditions in 
(O). □ 
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2.2 Choice of subtracted fields 

We now seek to find a choice for ip and X to use in Lemma 12.11 tliat leads to some optimal 
cancellation in the PDEs fl2.2l) - fl2.3l) . To this end, we first fix a gauge. Define 0o : ^ — to 
be the solution to (11. 9p . Note that 0o exists, is unique, and is smooth. We now fix a gauge 
in which = |log£:| (pQ. 

Lemma 2.2. It is possible to change gauges so that = |log£:| (pQ and so that B^{0) satisfies 
the Coulomb gauge, i.e. div B^{0) = and B^{0) ■ u = on dQ. 



Proof. The result is identical to that of Lemma 2.4 in [Ti]. □ 

With this choice of gauge, we now choose the ^p and X to work with. Define Hq : Q ^ 
to be the solution to (11.101) . Again, Hq exists, is unique, and is smooth. Define : ^ — ^ 
to be the solution to ffTTTD . Then we define Xq : Q ^ R"^ via ffrT2|) . which implies f lLTSj) . 
Finally, define ■j/'o ^ ^ ~^ 11^ to be the solution to (ll.Mp . This PDE is well-posed since (II. 3p . 
ffLQjl . and ffLTOj) imply that 

^.^ / aV0o-v^/io \ ^ r aV0o - v^hp r {hj-i)-u + i-v ^ f 

\ b J Jg^ b Jq^ b Jgn 

(2.5) 

The reason for defining ipQ and Xq in this manner is seen in the following lemma. 
Lemma 2.3. Let (po, Hq, Xq, and ipo be as defined in (II. 9p - (ll.Mp . Then 

aV0o - curl Xq = 6( V^/^o - ^o) • (2.6) 
Proof. Since divXo = 0, the definition of ipo (I1.14p implies that 

div(V^„ - Xo) = div , (2.7, 

SO that by Poincare's lemma. 



VV^o - ^0 = ^ h V X (2.^ 



for some x : — )■ M. Taking the curl of equation (12.80 yields 



Ax = -/^o-^TV^T-V(/)o + div(^ ) =0 (2.9) 



b • — ■ ^ ^ 

since Hq satisfies (ll.lOp . Taking the dot product of (12. 8p with the boundary normal yields 

- 9.x = V^X ■ ^ = ^ ■ ^ - ^^'^ ^ ^ • ^ = 0, (2.10) 

where r = i/"*" is the boundary tangent. This implies that x is a constant on the boundary, 
and then inside, and so the result follows since curlXo = Hq. □ 
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This property leads to a nice cancellation in the equations of Lemma 12. which allows 
to obtain an equation fl2.13l) of the form fll.51) (in a gauged version). 

Lemma 2.4. Let ipo, be as in flL14p and (IL12p and suppose {us, A^) are solutions to 
(11. ip in the = |loge| 0o gauge. Define = \\oge \ Z for 

Z:=VtPo-Xo (2.11) 

and also define 

/, = ^-|Z,|V/3|log£|Vo. (2.12) 
Then Ve = UeC-'^^"^'^^" /Vb and B, = - |log£| Xq solve 

(a + % |log£| P)dtV, = Ab^Vs + ^(1 - \Ve\^) + Vlogfe ■ VB.Ve + llZ, ■ W B^V^ + Vefe, (2.13) 

adtB, = V^h'^ + b{tv„VB,v,) + {\v,\^ - l)bZ, (2.14) 
in VL, along with the boundary conditions 

'h'=0 ondVL , , 

V_B^f£ ■ V = — ifgVlogft ■ u on dn. 

Here we have written h'^ = curli?^. 

Proof. We apply Lemma [2.11 with / = V'o and X = Xq. According to Lemma [2. 3 [ 

div(6(V^o - Xo)) = crA0o = a&0o- (2.16) 

Replacing $e by | log el 0o, this yields the cancellation of all the terms multiplying ive in the 
equation (12. 2p . which gives (I2.13p . A similar application of Lemma [2.31 gives (I2.14p . □ 



Remark 2.5. The vector field Z defined in Lemma 2^ is smooth and \\Z\\^2(^q\<oo. This 
follows immediately from its definition and the smoothness of the solutions to (ll.9p - (I1.14p . 



3 A priori estimates and compactness 

We introduce the pinned free-energy density 

g^iu, A) := ^ (^b + ^(1 - luff + Icurl Af^ , (3.1) 
and the free energy 

F,{u,A)= [ ge{u,A). (3.2) 

We are interested in proving a priori estimates and compactness on the magnetic field 
Bit. The choices of gauge and iPo^Xq give rise to nice properties for it. Some difficulty goes 
into controlling the terms in since the maximum principle does not hold for this mixed 
flow equation, and as such it is not guaranteed that \vs\ < 1. 
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Lemma 3.1. Let {vi;,B^) solve fl2.13p - fl2.15l) in the = |loge|0o gauge. Then the vector 
field satisfies 



adtdivBs = ab{iVe,dtVe) + a \\oge\ {\vef - l)60o + |loge| dt 



(3.3) 



in Q, and 

adtB, ■ u =\\oge\i\v,f-l) J -u 
on dQ. Consequently, for 1 < p < 2 and q = 2p/ (2 — p) 



(3.4) 



a\\diYBe{t)\\L, <a [ Vbdtv,{s) (Vb + Vb 

Jo ^ LI 



|1 - \Ve\"\\^, ) ds 



+ a\\oge\ / ||60o|Il^ ||(|^^£(-s)|^ - 1)||^2 
Jo 



and 



+ llogel (||1 - M)W, + ||1 - IvMfWJ (3-5) 



^PeW ■ '^IIlp(oq) < / \^oge\\\J ■ u\\^^(Q^)\\{\ve{s)f -l)\\,^^ff^^;^ds. (3.6) 



\LP{dn) 



Proof. The equation in (12. lip reads 



(3.7) 



Taking the dot product of this equation with the boundary normal v and applying the 
boundary conditions h'^ = Q and Vs^We ■ z/ = —^^v^Vlogh ■ v on f l2.15p yields 

adtB, ■ V = V^h', ■ V + h{%v,, ^ b,v, ■ v) + (\v,'\ - 1) |log£| (V^o - Xq) ■ u 

= {\v,f-l)\\oge\J-u, (3.8) 

which is (13. 4p . Taking the divergence of (13. 7p and employing Lemma [2. 3 j we find that 

adtdivB, + ab\\oge\(l)o = div {b{iv ^, V B,Ve) + \ve?bZ^). (3.9) 

On the other hand, taking (if^, ■) with f l2.13p yields the equality 

a{ive, dtVe) + /3 |loge| (w^, dtVe) = {iv^, Ab.v^ + V logfe ■ Vb^v^ + 2iZ, ■ Vb^v^) 

= diYiiVe, VB^Ve) + {iv,, V b^v^) ■ V logb + V {v^f ■ Z, (3.10) 

so that (again using Lemma 12. 3p 

/ 1 1 2 -|^\ 

ab{iVs,dtVe) + /3 \\oge\ bdt — ^— h ab \\oge\ {vel"^ (po = div {b{iv ^, V b^v^) + \ve\^bZ^). 



Equating ([S3]) and flXTT]) gives f lO]) . 



(3.11) 



Recall that by Lemma \2.2\ at time t = the vector field B^{0) satisfies div B^{0) = 
and -Be(O) ■ z/ = 0. The estimate fl3.6l) then follows directly from integrating f l3.4p in time 
from to t. For fl3.5p we first integrate fl3.3l) in time and then apply the Holder inequality 

< 



LP 



^2 ii<f iii<7 for q = 2p/{2 — p) to bound 
\\b{iVs,dtVe)\\^^ < \\b\v^\\dtv^\\\^p < \/bdtV^{s) ( Vh 



L9 



+ 



and 



60o(k(^)r-l)|L, <||&0oL,||(l^e(^)|'-l) 



Il2 



(3.12) 

(3.13) 
□ 



We will use this result with the following. Recall is defined in (13.20 . 
Lemma 3.2. For any (f£,i?e) (^noi necessarily solutions) it holds that for 2 < g < oo 



< 



^el IIL9 

for some universal constant C > 



F,{v,,B,) and ||1 - 11^,(9^) < Ce^''^^ F,{v,, B,) 



(3.14) 



Proof. Write p = By rewriting = ptu^ and remembering that < inf 6 < 1, we find 
that 



|Vp|V 



2\2 



2e^ 



2\2 



2^2 



<CF,{Ve,B, 



Notice that since p > we may bound (1 — p)^ < (1 — p^)^, so that 



\\l-ph.<Ce^/F,{ve,B,) and ||1 - p||^, < C/F^^^^. 
Interpolating between these bounds yields 

||1-pIIh^ < 111 - plli^' 111 - pII^i < Ce'-WFe{vs,B,) 
for any s G (0, 1). When s G (1/2, 1) trace theory then gives 

||1 - p||j^s-i/2(af^) < C ||1 - p||j:^.(f^) < Ce'^''^/F^{ve,B^). 



(3.15) 

(3.16) 
(3.17) 
(3.18) 



Then the bounds fl3.14p follow from fl3.17p and f l3.18p by using the embedding iF )■ L'^ for 
q = 2n/{n — 2r) first with r = s and n = 2 and then with r = s — 1/2 and n = 1. □ 

We will need the following inequality in order to derive some compactness results for B^. 

Proposition 3.3. For 1 < p <2 and 5 > sufficiently small there exists a constant C > 
so that 

\\A\\wyv-s,P < C (\\dwA\\^^ + ||curM||^p + \\A ■ ^Wlp^q^^ (3.19) 
for all A G W^^P{VL- M^). Moreover, 



LI 



< C (lldiv A||^, + ||curM||^, + \\A ■ v\\^,^q^^ 



(3.20) 



for q = 2p/{l + 5p) . 
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Proof. We employ the Hodge decomposition A = V0 + V'^ip for the solution to 

Ad) = divA inQ , , 

^ (3.21) 
Vcf) ■ u = A ■ u on dQ 

and ip the solution to 

Aih = curl A in Q , ^ 

^ (3.22) 
ip = on 50. 

The usual elliptic theory jLM] provides the existence such solutions satisfying the estimates 

<C\\cm\A\\^, (3.23) 



and 

lUI, 

\w 



< C (^||div All^j, + \\A ■ z/||^p(a^)) (3.24) 



for any S > sufficiently small. 
Then 

< + ||V0||,^i/p_a,p < ||V^||,4/2.p + (3.25) 

which together with ^23\f and (13:21) imphes flXTOD . The bound (13:201) follows from (13391) 
and the embedding W^^^'^'^ ^ L'^ for 

1 _ 1 1/p - 5 ^.^^ ^ ^ 2. (3.26) 



□ 



q p n 

We now turn to some a priori bounds on 1 — \vs\ and -B^. 
Proposition 3.4. Let (t>e,i?e) so/we (I2.13p - (l2.15p . Suppose that 

sup F,{vs,B,){s)+ [ [ab\dtVsf + a\dtB,f<K\loge\. (3.27) 
o<s<t Jo Jn 

Then the following hold. 

1. For any 2 < r < oo there exists a constant C depending on K such that 

sup ||1 - |t^£(s)|||^r < Ce^/'' ^J\\oge\ and sup ||1 - be(s)|||^,.(aQ) < Ce^''' \/\\og£\. 

(3.28) 



0<s<i 0<s<t 



2. For any 1 < p <2 and 5 > sufficiently small there exists a constant C > depending 
on K so that 

sup < C^\\ogs\{l + t\\oge\e^/'') (3.29) 

Q<s<t 

for q = 2p/{l + 6p). 
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3. For any 2 < q < 4 there exists a constant C > depending on K so that 

sup \\B,{s)\\^, < C^/\\oge\{l+t\\oge\e'/'^). (3.30) 

0<s<t 

4- There exists a constant C > depending on K so that 

[ \v,it)\'\Beit)\'<C\\oge\il + t'\\oge\'V^). (3.31) 

Proof. To begin, we note that Lemma 13.21 and fl3.27p imply that 

111 - Ivelhnn) < Ce'/''V\^\ and ||1 - \v,\\\Lridn) < Ce'^'' (3-32) 

for all r > 2, which yields (E2H1). 

To prove the second item, we let 1 < p < 2 and 5 > be sufficiently small. According to 
Proposition 13.31 we may bound 

\\Be\\wi/p-s,P < C (^\\div B^W^p + ||curl5^||^j, + \\B^ ■ i^\\LP(an)) ■ (3-33) 

We will estimate each term on the right hand side of this inequality using Lemma 13.11 with 
the bounds (13.271) and (I3.32p . The bounds (13.61) and (I3.32p together with Holder's inequahty 
provide the estimate 

a\\B,it) ■ u\\^,^g^^ < Ctlhgef^'/'^ (3.34) 
for q = 2p/(l + Sp). Similarly, (13. 5 p gives 

a||div5,(t)||^, <C 

<C\\ogef^e^/''{l + t). (3.35) 

Finally, the energy bound (I3.27P and Holder imply that for p < 2, 

1 1 curl 5, (t) 1 1 < C||curl5,(t)||^2 < ^v^k^. (3.36) 

We may then combine estimates (I3.34I) - (I3.36I) with (13.330 to deduce (I3.29P 

For the third item we now fix g G (2,4) and choose 1 < p < 2 and 5 > so that 
q = 2j9/(l + 6p) and 1/p — 6 > 0. Then (13.300 follows from (I3.29P and the embedding 

Now for the fourth item we utilize (I3.32p with r = 2g/(g — 2) > 2 to bound 

wivMrnmy < iiit'.(t)iiL.iii?.(t)iL, <cv^i^(i+tiiog£|£i/^). (3.37) 

Squaring this inequality and applying Cauchy's inequality on the right side then gives 

\\\ve{t)\ \B,{t)\\\% < C \\oge\ (1 + t^ llogel'e^/-?) < c \\ogs\ (1 + t^ \logef v^), (3.38) 
where in the last inequality we have used the fact that 2 < g < 4 implies 2/q > 1/2. This is 

(13311). □ 
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VbdtVe 
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We now parlay these bounds into convergence results. We begin by recalling a Lemma 
on compactness in space-time, due to Simon [Si] . 

Lemma 3.5. Suppose X, Y, Z are Banach spaces such that X CC Y ^ Z and 

II II <r r II 11^"^ II ll" 
IfIIy ^ IfIIx IfIIz 

for some 9 G (0, 1). Let 1 < pi,p2 < oo. Then each set bounded both in Lp^([0,T]; X) and in 
|yi,P2Q0,T];Z) IS pre-compact m Lp([0,T];Y) for all p < pi/{l - 9). 

With this lemma in hand we can deduce a pair of convergence results. 

Proposition 3.6. Let solve fl2.13p - fl2.15p . Suppose that 

sup F,{v,,B,){s)+ r I ah\dtv,\^ + a\dtB,f <K\\oge\ (3.39) 
o<s<T* Jo Jn 

for some fixed T* > 0. Fix 2 < r < 4. Then up to the extraction of a subsequence 

B. 



' ^ B, tn L\{0,n);L'{Q)). (3.40) 



Moreover, 



and 



V|log£ 
curl Bs 



curlE, zn L^(n x (0,T,)), (3.41) 



llog^l Jo Jn 



Proof. We will derive these convergence results by applying Lemma 13.51 As such we must 
first verify its hypotheses. We begin with the convergence of i?£ := y^|loge[. According 
to (Km of Proposition [33] 

sup ||5,(t)||^,/^_,,^ <a (3.43) 



0<t<T, 

On the other hand, the energy bound fl3.39l) implies that 





dtBM\%dt<C. (3.44) 



Hence the collection {B^}^ is uniformly bounded in 

L°°((0, T,); W^/P-^'^in)) n H\{0, n); L\n)). (3.45) 

According to the interpolation result of Theorem 4.3.1/1 of there exists a constant 
C > so that 

ll^elliy.,, < C \\Be\\L2 (3-46) 

for any 9 G (0, 1), where 

--6] ^nd-=^-^ + ^-. (3.47) 
p J q p 2 
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We chain this inequahty together with the embedding 

W'^iQ) ^ UiQ) for - = i - - with s > (3.48) 

r q 2 

to get 

\MLr<C\\B4'J,^^s.A\Be\\[. (3.49) 

for 

1 1-9 9 1-9/1 



,--6]. (3.50) 

By choosing l<p<2,S>0 small enough, and 9 G (0, 1) we can achieve any r G (2,4) in 
(I330|) . So, for any r G (2,4) we can apply Lemma [33] with X = W^/P-^'P{n), Y = L^(fi), 
and Z = L2(11) to get fl3^ . 



We now turn to the convergence of Q := curl _Be/A/|loge|. Because of fl2.14p we may 
write 

V C.=^5t5. r- 1 -{\Ve\ -1) r- 1 - (3.51) 

V|loge| v|loge| 

This equation, the energy bound f l3.39p . Holder's inequality, and f l3.28p then allow us to 
bound, for 1 < j9 < 2 

^* IIVCIIL <C j^^' {\\dtBe\\l. + ''^fjg^'p' + 1 1 < C(l + T,). (3.52) 
Similarly, the bound f l3.39p and Holder give, for 1 < p < 2, 



IICellL < C r < CT,. (3.53) 



Jo 



\loge\ 



We will estimate dtCs spatially in H-\n) = {H^{n)y. Let ip G L'^{{0,n); H^{n)). Then 
since (s = curl 5^ we may estimate 

I ' UdtCe= r I -V^^ ■dtB,<( [ ' PtBeWl.) ( [ * ll^llii ) , (3.54) 
which implies, with f l3.44p . upon taking the supremum over all such ip, that 

r WdtCeWl-. < r ptBeWl, < C. (3.55) 
Jo Jo 

Hence, from f l3.52p . fl3.53p . and (I3.55p we see that for any 1 < p < 2, the collection {(^} is 
uniformly bounded in 

L2((0, T,); W^'^m n H\{0, T,); H~\n)). (3.56) 

To apply Lemma [3.51 we fix 1 < p < 2 and apply Theorems 4.3.1/1 and 4.8.2 of to see 
that 

\\aw^.<C\\Ce\\wL\\Ce\\H-^ (3-57) 
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for 9 e (0,1), 

3 = 1-26, and - = + -. 3.58 

q p 2 ^ ' 

If we choose 9 = {2p — 2)/(3p — 2) G (0, 1/2), then s > and we may use the embedding 
f|g:ig]l with r = 2 to bound 

\\Ce\y<C\\Ce\V.\\Ce\\H-^- (3.59) 

Hence Lemma [33] is apphcable with X = W^^^iVL), Y = L'^iVt), and Z = H-\Q), and it 
provides the convergence (13.411) . The convergence of (13.421) is a consequence of (I3.4ip and 
the weak-L^(fi x (0, T^,)) compactness of dtB^/ y^jlogej that follows from the energy bound 

□ 



4 Energy analysis 

In this section we examine the evolution of the energy and show that it cannot increase too 
quickly in time. We first have the following formula for the evolution of the energy density 

Lemma 4.1. For any pair {u,A) (not necessarily solutions) it holds that 

dtQeiu, A) = bdiv{dtu, Vam) + curl (curl A (5t A)) 

- idtU,bAAU + ^{l-\tif)) - (dtA) ■ {V^ cml A + b{iu,VAu)). (4.1) 

Proof. The result follows from a direct calculation and the commutator identities 

VAdtU - dtVAU = iudtA (4.2) 

and 

{82 - iA2)idi - iAi)u - {di - iAi){d2 - iA2)u = tu curl A. (4.3) 

□ 

When we apply this to our solutions {vs,Be) in the $e = |log£:| (po gauge we are led (as 
seen in Section II. 4p to consider a modification of the free energy density given by 

~g,{u, A) := g,{u, A) + il^Mife/^, (4.4) 

where is as defined in Lemma 12.41 The evolution of the integral of this energy density 
involves a surface energy term, so we are led to define 

F,{u,A):= [ ~g,{u,A)+ [ l^tllvfe ■ j.. (4.5) 
Jn Jan ^ 
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Lemma 4.2. Let {v^iB;,) he solutions to fl2.13p - fl2.15l) in the = |loge|0o gauge. Then 

dtgeivs, B,) = div {b{dtVs, ^^B.Ve)) + cnilik'^dtB^) - ah \dtVef - cr \dtBe\^ 

+ hZ,-V{ve,B,). (4.6) 

In particular, this implies that 



dtF,{v,,B,)+ I ah\dtV,\' + a\dtBf= I hZ, ■ V{v,, B,). (4.7) 
Jn Jn 

Proof. We plug the equations of Lemma 12.41 into Lemma 14.11 to find 

dtgeive, B^) = hdiv{dtVs, Vb,v,) + cnT\{cni\ B,{dtB^)) 

- {dtVs, {a + i \loge\ P)hdtVs - Vh ■ Vb^v^ - 2thZ, ■ Vs^Ve) 

- (dtB,) ■ {adtB, - {\vef - l)hZ,) + {dtv,,v,)hf,. (4.8) 

The equation (14. 6 p follows from an expansion of these terms and the equality (see Lemma 
2.12 of [E]) 

{dtVe,2iy BVe) = V{Ve,B,) + (1 - \v,\^)dtBe. (4.9) 

The equation (14. 7p follows from integrating (14.60 over VL and applying the divergence theorem 
and the boundary conditions fl2.15p to get 



/ div {h{dtv,, Wb.v,)) + cm\{h'^dtB,) = [ hidtv,, ~v,\/ hgh ■ u) 
Jn Jan ^ 



/ lt^i^-^V6 ■ i.. (4.10) 
ho. 4 



□ 



Since our modified energy is not positive definite, we must prove a result that controls 
in terms of and a negligible error. This is the content of the next lemma. 

Lemma 4.3. For e sufficiently small we have the estimate 

F,iv,,B,) < {l + Ce'/^)F,{ve,B,) + CE'/'' (4.11) 

for a constant C > 0. 

Proof. Write x '■= V6 ■ u. We apply Cauchy's inequality to find 



hi\Vs\'-l)fe 



|2 



<|F,(t;„5,) + | /" \f,\' <'-F,{ve,B,) + Ce\\ogef. (4.12) 
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For the boundary term we estimate with Holder 



(kl'-i)x 



dn 



By Lemma [3.21 we know that 



\vA - II 



Li{dn) 



< Ce'/'^./F,{v,,B,) 



for any 2 < g < oo. We may then rewrite 



l-\v,\^ = -{l-\v,\f + 2{l-\v. 



(4.13) 



(4.14) 



(4.15) 



in order to bound 



\v^\ -1 



<C \\\v. 



\vA - II 



<Ce'/\F,{v„B,) + y/F,{vs,B,)). (4.16) 
We then combine f l4.13p and f l4.16p and again use Cauchy to get the estimate 

(kl'-i)x 



'an 4 

Now from (l4T2il and (14T71) we know that 



<Ce'/\FM,Be) + l). 



(4.17) 



F,(^;„B,) = F,(t;„5,) + 



Ki - b.l ) 



(Kl'-i)x 



an 



> F,iv„B,) - -F,iv,,B,) - Ce \\oge\' - Ce^/\F,{v,,B,) + 1) 

> (1 - Ce^/^)F,{ve, B,) - Ce^/^ (4.18) 

when e is sufficiently small. This yields the desired inequality. 

□ 

As we mentioned, the control of the energy growth comes from the product estimate from 
|SSlj . which we now state in our context. 

Proposition 4.4. Let (v^^Bi;) solve fl2.13p - fl2.15p and suppose that 



sup Fe{Ve,Be){s) + 
0<t<T, 



ab\dtVef + a\dtB,f < C|log£| 



(4.19) 



'0 Jn 

for some fixed T* > 0. Then the following hold, up to extraction of a subsequence. 



1. There exist /i G L°^([0, T,]; A^(fi)) and V G L'^{[Q,T^]]{M{VL)f) such that dtH + 
curly = and I2{v„ B,) fi, V{ue,Be) ^ V as e ^ in {C°'\n x [0,T,]))*. 
Here we have written = (C°(r2))* for the space of bounded Radon measures. 
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2. For any [^1,^2] ^ [0,r,], any Y eC^i^x [^1,^2];^^), G C°(^] x [ti,t2]), we have the 
bound 



t2 



ti Jn 



< liminf 



i / / / ,,|2 



£^0 I log e I ifi ^ Jn 



1/2 



(4.20) 



c?. r/ie mapping t 1— )■ {fi{t),^) is in if^([0,T^,]) /or any ,^ G C^(ri), anc? m particular 



Mt2),0 - (M^i),OI < Cv/tJ^liminf ^_ 

V|loge| Vti 



i2 



\dtVe 



1/2 



(4.21) 



/or any [^1,^2] ^ [0,T]. 
Proof. To begin we note that 



|Vt;e| < IVb.VsI + \Vs\ Ws 



(4.22) 



so that f l3.3ip of Proposition 13.41 and f l4.19p . combined with the fact that b is bounded below 
by a constant independent of e, imply that 



2n2 



o<t<T. 2 



4e2 



Jn 



\dtVef < C \\oge\ 



(4.23) 



This means that Theorem 3 of |SSlj is applicable and provides the results of item 1 for 
V^(f£,0) and /i(fe,0). However, the definitions of ( fll.191) and fll.201) ) imply that 



and 



^^{v,, B,) - fi{v,, 0) = curl ((1 - \v,\^)B,) 



V{v„B,) - V{v„0) = dt ((1 - \vs\')B,) . 



(4.24) 



(4.25) 



The bounds on and on B^ f l3.30p imply that V{ve,Bs) — ^(^£,0) — )■ and fi{vs,0) — 
fi{ve,Bs) — )■ in the dual of 1^^'°° and hence in (C°'^)*, so the results in item 1 hold as 
stated. 

For item 2 we note that Theorem 3 of |SSlj also yields the estimate 



ti Jn 



V ■i)Y 



< 



1/2 / ^ \ 1/2 



(4.26) 



'nx[ii,t2] / \Jnx\tiM] 
where vy and are the defect measures of L'^iVL x [0,T*]) convergence of 

and 



Vllogel \/\^oge\ 



(4.27) 



respectively. However, 



\^B^Ve ■ Y\^ = \Vv, ■ Y\^ - 2{lVe, Vv,) ■ Y{B, ■ Y) + \v,\^ \Be ■ Y\ 



(4.28) 
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so the convergence results in Proposition 13.61 and the bounds fl3.28p of Proposition 13.41 guar- 
antee that the defect measure of L'^{Q x [0,T*]) convergence of 



coincides with z/y. Then (14.201) follows from f l4.26p since 



(4.29) 



1/2 / /■ \ 1/2 



nx[ti,t2] / \Jnx[ti,t2] 



<liminf-^f r / \VB^v,-Yf T / \ijdtvA ' . (4.30) 
e->o |loge| VAi Jt^ Jn J 

Item 3 follows easily from item 2 as in Theorem 3 of |SS1] . □ 
We then deduce the control on the energy growth, following the same method as in [Ti] . 

Theorem 4.5. Let {vs, -Be) be solutions to ( ]2.13p -( !2.15p in the $e = |loge| 0o gauge. Suppose 
that the initial data (ve(0), Be{0)) are well-prepared in the sense of (11.220 . Fix Cq > 0. Then 
there exists a constant Tq = To(Co) > such that, as e — )■ 0, 

F,{ve,Be)it) < Feiv„Be)iO) + Co\\og6\ for all t E [0,To] (4.31) 

and ^ 

r I ah \dtv,\^ + a \dtB,\^ < Co \hge\ . (4.32) 
Jo Jn 

Proof. The proof is essentially the same as Theorem 2.16 of [Ti], so we will present only a 
sketch of the idea. For full details see [Ti]. For any t > 0, consider the two conditions 

F,{ve,B,){s) < F,{ve,B,){0) + Co \\oge\ for all s G [0,t] (4.33) 

and ^ 

[ [ ab \dtv,f + a \dtB,f < Co \\oge\ . (4.34) 
JO Jn 

Define 

% ■■= sup{t > I conditions (SSS]) and (S31 hold}. (4.35) 

The smoothness of [ve, B^) and the well-preparedness of the initial data guarantee the exis- 
tence of a time t^ > (depending on e and Co) such that both conditions hold for t^. Hence 
7e > for each e. We show that actually 7^ > Tq for some Tq > as e — )■ 0, thereby proving 
the theorem. Suppose, by way of contradiction, that 

liminf7^ = 0. (4.36) 

e-5-O 

We may suppose, up to extraction of a subsequence, that 7^ — )■ as £ — > 0. 
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Rescale in time at scale 7e by defining We{x,t) = Vg^Xj'j^t) and Cs{x,t) = Bir{x,'yet)- By 
the definition of 7^, the inequalities 

F,{w,, C,){t) < F,{we, C,)(0) + Co \\oge\ for all t e [0, 1] (4.37) 

and ^ 

- / j ah \dtw,\^ + a \dtC,\' < Co \\oge\ . (4.38) 
% Jo Jn 

both hold, but at time t = 1 one of the inequalities must be an equality since We and Ce are 
smooth. However, these inequalities and Lemma [4.31 provide a bound for Fe(i(7£, C^), which 
in turn allows us to apply Proposition 14.41 to deduce, since 7^ — > 0, that V{ws, Cs) — )■ and 
n{we,Be) — )■ A*(0). This information and Lemma [4. 2[ rescaled in time at scale 7^, then show 
that neither inequality fl4.37p nor f l4.38p could be an equality at t = 1, a contradition. Indeed 
from V^Wir, Cg) — )■ and (14. 7p properly scaled in time we deduce 

F,{w„C,){l) < F,{w„C,m + o{\\oge\), (4.39) 

so there cannot be equality in (I4.37p . On the other hand since the limiting vortices do not 
move (from \^ = 0) by the F-convergence lower bound (11.230 we must have 

F,{we,a){l) >7lJ2Ka^i^)) \log6\ + o{\log6\) = F,(w„ C,) (0) + o( |log £| ) . (4.40) 
Combining with (14. 7p again we deduce 

- r / ab\dtwe\' + a\dtCe\' = o(|log£|), (4.41) 
7e Jo Jn 

and hence (I4.38P cannot be an equality either. □ 

We now improve the comparison between and F^ from the result of Theorem 14.51 

Corollary 4.6. Suppose the hypotheses of Theorem \4.5\ Then as e ^ 0, 

F,{ve,Be){t) < F,{ve,B,){0) + Co\loge\+o{e'/') for all t e [0,To]. (4.42) 

Proof. This is a direct consequence of Lemma 14.31 the well-preparedness and the bound on 
F^ provided by Theorem 14.51 

□ 



5 Convergence results 

From here on we assume the various assumptions made in the introduction, in particular 
(ll.22p . so that there are n initial vortices of degrees ±1. We now fix Cq > so that 

Co < TT ( inf 6(x) ) (5.1) 

and apply Theorem 14.51 to get To > so that the conclusions of the theorem hold on the 
interval [0,To]. Then Corollary 14.61 provides an estimate for Fs{ve,B^) for all t E [0,To]. 



29 



5.1 Vortex trajectories 

Our first result allows us to define the n vortex paths in the time interval [0,To]. 
Lemma 5.1. Let 

70 = min{|ai(0) - aj(0)| \ i ^ j} U {dist{ai{0) , dQ)} . (5.2) 
Fix < 7* < 7o . Then there exists aT^, = T*(7*, Co) with T* G (0, To] so that 

1. The space-time Jacohian {fi{vs, Be),V{ve, B,)) ^ {fi,V) m {C^^^{VL x [0,To]))*. 

2. There exist functions ai G H^{\Q, T^\,VL) , i = 1, . . . so that for all t G [0, T,,] 

mm{\ai{t) - aj{t)\ | i ^ j} U {dist(ai(t), Sr^)} > 7*, (5.3) 

n n 

fi{t) = 2nJ2dMKit), and V{t) = 27r ^ ci,(0)d,^(t)5,,(t). (5.4) 

i=l i=l 

Proof. The bounds provided by Theorem 14.51 and Corollary 14.61 allow us to apply Proposition 
14.41 to deduce the first item. The second item would be standard if 6 = 1 since then the 
quantization of the energy would prevent the nucleation of new vortices. However, when 
b ^ 1 the energy is not quantized and so we must look for another structure to prove the 
result. 

We first recall that by the third item of Proposition 14.41 the mapping t 1— )■ {il!,^{t)) is 
(and hence continuous) for any fixed ip G C^{Q). Fix < 6 < 7r(inf 6) — Cq. By the 
smoothness of b, there exists 77 > so that 

\b{x) - &(a^(0))| < 5/{n7i) for all x G U^=iS(ai(0), r^). (5.5) 

Standard "Jacobian estimates" (see |JS2] or Chapter 6 of jSS3] ) combined with the energy 
upper bound and the lower bound on b, give the quantized structure of the measure /i(t), 
i.e. 

n{t) 

fi{t)=27fJ2MtKm (5.6) 

i=l 

for |(ij(t)| > 1. Moreover, the F— convergence lower bound fll.23p . when combined with 
f02]) . shows that 

n(t) n 

7r^|ci,(t)|6(a,(t)) <7r^6(a,(0)) + Co. (5.7) 

i=l 1=1 

Let 7i = min{(7o — 7*)/2,?7}. For i = 1, . . . ,n let ipi ^ be a function supported 

in S(aj(0),7i) so that ipi{x) = 1 for x G -B(aj(0), 7i/2). Since t 1— )■ {tpi, ii{t)) is continuous 
and {ilji,^{0)) = 2ndi{0), we see that for t <T^ with T, sufficiently small at least one of the 
points aj(t), j = 1, . . . ,n(t) must be contained in the ball -B(aj(0), 71/2). Up to relabeling 
the indices, we may assume that aj(t) G i?(aj(0), 71/2) for i = 1, . . . ,n. This implies that 
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n{t) > n ioi t < T^. Now suppose by way of contradiction that n{t) > n for some t < T^^. 
Since aj(t) G 5(0,(0), r^) for i = 1, . . . , n, the bound (I5.5p imphes that 



TT 



J2b{a^{t))-b{am) 



1=1 



< 5. (5. 



Plugging this into ( \5.7\i . we deduce the bound 



n{t) 

Co + S>7i \d,{t)\ b{ai{t)) > {n{t) - n)7r(inf 6). (5.9) 

i=n+l 

The choice of S and Co then imphes that n{t) — n < 1, hence n{t) = n for Q <t <T^. The 
structure of V follows from the structure of n and the equation dtii + curl V = 0. 

Returning to the continuity of t h- >■ {ipi, fJ'it)), we find that di(t) = di{0) for i = 1, . . . ,n. 
Since aj(t) G -B(aj(0), 7i/2), we know that 

mm{\ai{t) - aj{t)\ | i ^ j} U {dist(ai(t), 9^)} > 7,. (5.10) 

Finally, to see that G if ^([0, T^^]; fi) we use the fact that the mapping t {ip,fi{t)) is in 
for the functions 

ip^x) = {x ■ e)ipi{x),i = 1, . . . , n (5-11) 

for any unit vector e G M^. □ 

Remark 5.2. Since the degree of the i*^ vortex does not change for any time in [0,T^], we 
may consolidate notation and write only di in place of di{0). 



5.2 Normalized energy density 

We can now show that up to the extraction of a single subsequence, the modified energy 
density converges for all t G [0,T^,]. 



Proposition 5.3. There exists a subsequence so that 

U'"e,B,){t) 



\\oge\ 



v{t) (5.12) 



weakly-* in {C^{VL))* for all t G [0,T*]. For each t G [0,T*], v{t) is a measure, and v G 
L°°([0, T^,]; ((])), where hAi^l) is the space of Radon measures on Q. Finally, 

^J^^Il^ ^ y(t)dt (5.13) 

|loge| 

weakly in the sense of measures on Q x [0,T*]. 

Proof. The proof of f l5.12p is the same as Proposition 4.4 of [Ti]. The fact that z/(t) is a 
measure, the inclusion u G L°°([0, T*]; A^(f2)), and the weak convergence (15.131) all follow 
easily from the energy bound of Corollary 14.61 and weak compactness of measures in f2. □ 
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Since we do not have precise control of the energy evolution yet, we cannot derive the 
exact structure of the measure v. We can however, provide a lower bound that follows 
directly from fll.23p . 

Lemma 5.4. It holds that 

v{t)>^-W)\- (5.14) 

5.3 Stress-energy tensor, etc 

We define the stress-energy tensor associated to (f^, B^) by 

fh , ,2 . , ,2x9 1 , 

Te = bVB^Ve ® VB^Ve - I2x2 l"^ B,V,\' + —(1 - l^;,!')^ - - |curl5, 

+ /2X2|^^C^(&A)1- (5-15) 

Note that 

T, = hWB^Ve ® VB^Ve - h>,2 {9e{v„ B,) - |curl5,|') . (5.16) 

The next result provides for the convergence of the (normalized) stress-energy tensor 
along with a couple other quantities. 

Lemma 5.5. Let he given by f IS.lSp . Then up to the extraction of a subsequence, as £ — )■ 
we have 

T, (5.17) 



|loge 

b{dtVe, Vbv, 
llogel 



^p, (5.18) 



and 



h\dtVe\^ \dtBe\^ 

— r + ^Ti — (5-^9) 

|log£:| |ioge| 



in the weak sense of measures on Q x [0,T^,]. 

Proof. All the terms defining are seen to be controlled by the energy density Therefore, 
Cauchy's inequality and the energy bounds of Theorem 14.51 and Corollary 14.61 show that 

^* ' ITJ b(dtVe,VB,v,) bldtVef IdtBf ^ 

\ s\ ^ V t e, — £^ ^ +^ ,/ I < C. (5.20) 



Jn 



\\oge\ |loge| \^oge\ |log£ 



Compactness in the sense of measures follows directly from these bounds. □ 

If we knew more information about the evolution of the modified energy F^, we could say 
more about the structure of the limiting measures T, p, and (. Instead we will resort to a 
Lebesgue decomposition of all the measures with respect to the vorticity measure fi, which 
we recall is equal to 27r diSai{t)dt. 
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Lemma 5.6. We can decompose 



V = uo{t)dt + ^ z/i(t)5„,(t)rft, T = To + ^ Ti{t)6a,(t)dt, {5.21] 

i=l i=l 



n 



C = Co + ^Q{t)5a,{t)dt, andp = po + ^pi{t)5ai(t)dt, (5.22) 

i=l i=l 

where Tq, (q, and \po\ are mutually singular with /i as measures on Q x [0,T*] and z/o(t) is 
mutually singular with fi{t) as measures on Q for a.e. t G [0,T*]. We have 

z/oGL°^([0,T,];A1(n)), (5.23) 

where A^(fi) is the space of Radon measures on Q. For i = 1, . . . ,n we also have that 

z/, gL°°([0,T,];M). (5.24) 

Finally, we have that i^oXoXiit) ^ 0, and the quantities Uiit) obey the hounds irb^aiit)) < 
< Trn ||6||j;^oo + Co for i = 1, . . . ,n, where Cq > is the constant chosen at the beginning 
of Section O 

Proof. We perform a Lebesgue decomposition of the measures u, T, (, and p with respect 
to /i in order to write 

u = uo + Uf,, T = Tq + T^, C = Co + Cm> and p = Po+ p^, (5.25) 



At! 



where z/o, |To| , Co, and |po| are mutually singular with as measures on x [0,T,] and 
Tfi, C^t; and p^ are absolutely continuous with respect to fi. According to Lemma ISTT] we may 
write fi = 2TrJ2i=idiSa^(t)dt, so that fi is supported along the curves {{ai(t),t) \ t G [0,^,,]}, 
which are parameterized by functions. This allows us to further decompose 



iyi{t)5a,(t)dt, T^ = Y1 Tr{t)5a,it)dt, (5.26) 

i=l i=l 

and 

n n 

Cfi = YCi{t)Sa,{t)dt, Pf, = Ypi{t)5a,(t)dt, (5.27) 

i=l i=l 

where z/j, T,, Ci, Pi are the Radon-Nikodym derivatives of z/^, T^, Caj? P^t along the curve 
{(ai(t),t) I t G [0,T]}. The decompositions and imply (IF:^ and the T decom- 

position in f l5.2ip . 

To finish the z/ decomposition in fl5.2ip we must extract more structure from z/o. Proposi- 
tion |53] implies V = i'{t)dt, where the mapping t h-> i>{t) is in L°°{[0,T^]] Using this 
and the above decomposition, we know that z^o G L°°([0, T,,]; A^(f2)) and z/j G L°°([0, T,,]; M) 
for i = 1, . . . ,n. This implies the decomposition of i> in fl5.2ip as well as the i> inclusions in 
and (^M>- 

The bounds i^o{t) , (o(t) , Qit) ^ are trivial. The lower bounds z/j(t) > 7ib{ai(t)) for 
i = 1, . . . ,n follow from Lemma [5^ The upper bounds z/j(t) < nn \\b\\^cx> + Co for i = 1, . . . ,n 
follow from the energy upper bounds of Corollary 14.61 and the well-preparedness assumption 

(Ol. □ 
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We now prove some estimates of |po|, l^ol, and \a cwc\ B^:dtB^ \ in terms of Co and uq. 
Lemma 5.7. It holds that 

\To\ < Ciyo{t)dt. (5.28) 
Let T] G C°([0, T^,]; M) satisfy rj{t) > for all t G [0,T*]. Then we may estimate 

I I , Co , VQ{t)dt 
2a 7] 

as well as 

\acm\B^dtB^\<ri— A . (5.30) 

' ' 2 ?7 

Proof. From (I5.16P we see that \T^\ / \loge\ < Ccjeive, B^) / |loge|, and passing to the hmit 
reveals that |T| < Cv. The decomposition of T provided by Lemma 15.61 implies that 
|T| = |To| + Yll=i \Ti{t)\5a^(t)dt, with |To| mutually singular with /x, and hence also with 
^r=i ^a,(t)dt and Yll=i ^i{i)^ai{t)dt. The estimate fl5.28p then follows. 

Now let T) G C°([0, T^,]; M) satisfy rj{t) > 0. The Cauchy-Schwarz inequality allows us to 
bound 

|log£:| ~ 2 |log£:| rj 2|log£:| 
Passing to the limit, we find that 

bl < + -• (5.32) 
la rj 

The bound (I5.29p follows from this, the decompositions of C and v provided by Lemma 15. 6[ 
and the fact that \p\ = \po\ + J2^=i \Pi\^ai with |po| mutually singular with J2^=iCi^a, and 

Similarly, we may bound 

\acmlBAB^\ ^ 7] a\dtB,f ^ a |curlE,|^ 



|log£:| 2 |log£:| rj 2|log£ 

which implies, upon passing to the limit, that as measures 



\acm\B,dtB^\ < ^ + — . (5.34) 
' ' 2 rj 



Then fl5.30p follows from the decompositions of Lemma 15.61 and (I3.42p of Proposition 
which implies that \a cml B^:dtB^\ is mutually singular with /x. □ 

We now compute the divergence of and relate its hmit as e — to the Lebesgue 
decomposition of p, z/, T. 

Lemma 5.8. The stress- energy tensor satisfies 

divT, = ab{dtV„VB,v,) -a curl B.dtB^ - ^^^^^V{ve, B,) - fi{v„ B,)bZ^ 

- iuve. B,) - V log 6 + [/3 |iog£| hdtB, + hWfe - fsVb] . (5.35) 
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(5.36) 



Dividing by |loge| and passing to the limit e — )■ 0, we have that 

div T = ap- - fibZ^ - (^u - V log 6 - a curl B.dtB^, 

where is the vector field given by Proposition \3.(A We have the decomposition 

n 

dwT = So + J2S^it)^Mt)dt^ (5-37) 



i=l 

where 

triTo) \ 



So = apo - — j V log 6 - cr curl B^ dtB^ (5.38) 

with \So\ mutually singular with fj,, and 

( tr(Ti) \ 

Si = api - l3Tcdib{ai)dj- - 2Tcdib{ai)Z^{ai) - f z/j ^ j Vlog6(ai). (5.39) 

Proof. A direct calculation reveals that 



divT, = {bAB_.v, + b\6-\l - \vf ) + Vb^v, ■ Vb, Vb^v,) - KAKiv,, Vb^v^) + V^h'^ 

4£2 



(5.40) 

By plugging in the equations fl2.13p - fl2.14p and using (14. 9 p we may rewrite 

{bAB^Ve + b\e-^il - \v,f) + V B^Ve ■ Vb, Vb^Vs) - KiKiVe, ^B^Ve) + V^K)^ 

= ab{dtv,, Vb^v,) - /^I^Qg^l V (t;„ B,) - aKdtB^ - fxiv,, B,)bZ^ 

-.,.o.,^^a-,.,Vv(<^)..,, 

where we have employed the identities 

{-2zVB,Ve ■ bZ„ Vb^v,) = 2{d^^ve, id^''Ve)bZj (5.42) 

and 

d^;,!' - l)h', + 2{df^Ve, id^^V,) = - CUTl{zVe, Vb^v,) - h'^ = -fl{Ve, B,) (5.43) 

to identify the —fi{vs, Bs)bZ^ term in (15.411) . From (15.150 we may calculate 

tr(T,) = ^ |2 ^ (|^^|2 _ ^^^^^^ 
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which imphes 



|2 



tr(T, 



Vlog6 -W^e.B,) + ^-^ - {\v,\' - 1)/,V6. (5.45) 



Combining fl5.40p - fl5.4ip with fl5.45p then yields (15. 35 p . Equation f l5.36p follows by dividing 
f l5.35p by I log £ I and passing to the limit e — )■ 0, using the fact that |log£:|^ (1 — \ve\^) — )■ 
(by the upper bound on the energy) as well as Lemma 15.51 and Proposition 13.61 to identify 
the structure of the limit. The decomposition fl5.37p with 5*0 given by fl5.38p and Si given 
by (15.391) follows by decomposing the terms on the right side of fl5.36p according to Lemmas 
15.11 and |5.6[ noting that (13.421) of Proposition 13.61 implies that cr curl iJ^.^t-B;;;'- is absolutely 
continuous with respect to Lebesgue measure on i7 x [0,T^,], and hence mutually singular 
with ji. □ 

Now we show that l^ol can be controlled in terms of VQ{t)dt and Co- 
Lemma 5.9. Let Sq he as in Lemma \5. <g[ For any rj G C'^([0, T^,]; M) satisfying 7]{t) > 0, we 
have the estimate 

\So\ < vCo +(c + uo{t)dt. (5.46) 

Proof. Using the definition of So from Lemma 15.81 along with estimates of Lemma 15. 7[ we 
may bound 

l^ol < a|Po| + iiyo+\To\)Vhgb+\acuT\B,dtB^\ 

+ + Cuo + — + = vCo+ [C+ z/o, (5.47) 

which is fl5.46p . Here we have used the fact that b is smooth and bounded below, which is 
guaranteed by (11. 2p . □ 



6 Dynamics 

We now use the convergence results of the last section to derive the dynamics of the vortices. 
We begin with a result that allows us to relate pi to the vortex velocity, hi. This result does 
not constitute the full dynamical law for since we do not yet know the value of Pi or z/j. 

Proposition 6.1. For i = 1, . . . ,n, it holds that 

p^{t) = -Vi{t)di{t) (6.1) 

for a.e. t E [0, T^]. 
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Proof. Fix < 77 < 7* (with 7^, given by Lemma [5.11) and a smooth vector field Y : [0, T*] — ?■ 
M?. Let be a smooth function with support in _B(0, 1) x [0, T*] so that V0(O, t) = Y(t) for 
all if: G [0,T^,]. Fix i G {1, . . . ,n} and let a] : [0,T^] — )■ be a smooth mollification of the 
path ai{t), with 7 > the mollification parameter. Define ip^{x,t) = ri(f){{x — a'J(t))/r],t) 
and ip{x, t) = ri(f){{x — ai{t))/'rj, t). 

Then the energy evolution equation f l4.6p implies that 



^,(T,) f (K(r.)i^-i) ^^ /• ^,(0) /• (K(o)i^-i) ^^ 

^^77; r+ / TTi < — yo-u- I I Ip^ — — . — Vo-z/ 

llog^l Jdn 4|loge| |log£| Jq^ 4|loge| 

in 



|log£| |loge| I Jo \\oge\ 



[ * / ^,V{v,,B,)-bZ+ [ * [ dt^,-^. (6.2) 
Jo Jn Jo Jn |loge| 



We may pass to the limit e — )■ in the last equation by using Propositions 13. 4[ 13. 6[ and 15.31 
along with Lemmas 15.11 15.51 to get 



/ ^^z/(T,) - [ ^^u{0) + [ [ ^■yC= [ [ -VV^^ ■ ip + cmlB,dtB^) 
Jn Ju Jo Jn Jo Jn 

+ 11 ^^V-Z+ I * / dtij^u. (6.3) 
Jo Jn Jo Jn 

Now, dtip-y{x, t) = rjdt(t){{x — a]{t))/ri, t) — a] ■ V0((x — a]{t))/ri, t). Letting 7 — )■ 0, using the 
boundedness of all of the measures involved given by Lemma 15.61 and employing dominated 
convergence, we deduce that 

[ MT*)~ I M0)+ r I i^C= r [ -Vij-{p + cni\B,dtB^) 
Jn Jn Jo Jn Jo Jn 

+ [ * [ ^V-Z + [ * [ {T^dMx - a,it))/r],t) - d,{t) - Vcjyiix - a,{t))/r],t))iy. (6.4) 
Jo Jn Jo Jn 

Note that here cij is the time derivative of aj G H^{[0,T^]) hence is an L^([0, T^,]; M^) func- 
tion, while z/ is L°°([0, T^,], A^(f2)). Hence the product ctji/ makes sense as an element of 
L2([0,r,];R2®A^(fi)). 

Recall that HV'IIloo + llf^tV^llLoo ^ Ci] and that the support of ip lies in a r;— neighborhood 
of the path aj(t). Passing to the limit — )■ reveals that 

iju{T,) - [ ^z/(0) + [ [ i'C^O (6.5) 



Jn Jo Jn 

and ^ 

[ i^V ■ Z + ridt(j)i{x - ai{t))/r],t)i^ ^ 0. (6.6) 
/o Jn 

Since Vip is supported in B{0,ri) and HV-i/'H^oc < C < 00, we also have that 

f ' f _ ■ curl B, dtBi 0. (6.7) 
Jo Jn 

37 



We may then pass to the hmit — i- in (16. 4 p and utihze Lemma 15. 6[ and the fact that 
V0(O,t) = V?A(ai(t),t) = Y{t) to deduce that 

= / Y-ipi + um). (6.8) 
Jo 

This resuh holds for any choice of F G C°^([0, T^]; M^), which imphes that 

Pi = -Vihi (6.9) 
for a.e. t G [0,T^,]. This is the desired result. 

□ 

Now we can use a similar argument to show that Ti{t) = for a.e. t G [0,T*]. 

Proposition 6.2. For a.e. t G [0,T=|,] and i = 1, . . . ,n it holds that Ti(t) = 0. 

Proof. Let ip G C^(M^;R) be such that supp('?/') G ^(0,7=,,), where 7^, is given by Lemma 
EH and so that ^{x) = 1 on 5(0,7,/2). Let K, G C0([0, T,]; M2x2) for i = 1, . . . , n. Define 
the vector field 

EM =,±m) . t^^^^ ( --^^) . (6.10) 

Since Lemma O says that G if ^ ^ C^'^/^, we have that E G C°(r] x [0,T*];M2) and 
DE G C°(fi X [0,T=|,]; M^^^). Note that, unlike in Proposition 16.11 we do not need continuity 
of dtE, so we do not have to use a smoothing of the vortex paths. 
According to the divergence theorem, we have 

' ■.DE= / j—^-E. (6.11) 



Jn 



|loge| Jo Jn |log£ 



By the above continuity results, we may pass to the limit in (16. lip and employ the decom- 
positions of T and div T, given respectively by Lemmas 15.61 and 15. 8[ to see that 



- [ * [ To -.DE - [ * VT,(t) :DS(a,(t),t) 
^0 Jn Jo 

= r [ So-E+ r*f^5.(t)-S(a.(t),t). (6.12) 
Jo Jn Jo 

By construction S(aj(t),t) = and DE{ai{t),t) = Ki{t), so (I6.12p becomes 

/ J2T,{t) : K,{t) = - To-.DE + So-E. (6.13) 

Jo Jo Jn 

The vector field E satisfies ||S|| roofovfo r n — ^V- This implies that 



L°°(nx[0,T.]) 
5*0 ■ ^ 



Jn 



<Cr] I / ISol ^ as ^ (6.14) 
Jn 
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since |S'o| has finite mass. Also, since Tq is singular with respect to fi and DE satisfies 
||DS||^oo(Qx[o T.]) — ^ supp(DS(-, t)) C U"^^-B(aj(t), 77), we have that 

' To : DE ^ as ri ^ 0. (6.15) 



'0 Jn 

Hence, taking the limit 77 —> in fl6.13p . we find that 

/ VT,(t):ir,(t) = 0. (6.16) 
-^0 i=i 

Since the Ki were arbitrary, we immediately deduce that Ti{t) = for a.e. t G [0,T^,]. □ 

Now we can deduce some estimates for Si, as defined in Lemma [5.81 

Lemma 6.3. Let Si, i = 1, ... ,n be as in Lemmas \5.S\ Then Si G L^([0, T^,]; M^) for each 
i = 1, . . . ,n. Moreover, if Yi G L^{[0, T^,]; M^) for each i = 1, . . . ,n, then for any t G [0, T*] 
we have the estimate 

for some C > 0. 

Proof. The result of Proposition 16.11 when combined with the definition of Si given by (I5.39P 
in Lemma 15.81 and the vanishing of Tj given by Proposition 16.21 implies that 

Si = —avihi — /3ndib{ai)dj' — 27r dib{ai)Z^{ai) — ViV \ogb{ai). (6.18) 

Lemma 15161 implies that Ui G L°°{[0,T^:]), Remark [231 shows that Z is bounded, and assump- 
tion (11. 2p provides the boundedness of b and V log b; then since we know from Lemma 15.11 
that di G L2([0,r,];R2), we find that Si G ^^([0, T,]; R^). 

We now turn to the proof of (I6.17p . assuming initially that Yi G C°([0, T*]; M^) for each 
i = l,...,n. Let ip G C^(R^;M) be such that supp('?/') G 5(0,7*), where 7* is given by 
Lemma [5. 11 and so that tplx) = 1 on i?(0,7*/2). Define the vector field 

n 

S(x,t) = 5^y,(t)V^(x-a,(t)). (6.19) 

i=l 

For each fixed t we have supp(H(-, t)) C U"^^i?(aj(t), 7*), and the choice of 7* implies that 
aj{t) ^ B{ai{t),'y*) for i ^ j. Since G C°'-^/^ by Lemma [5H we know that S G C^{Q x 
[0,r*];R2) and DE G C°((] x [0, T*]; R^^^). if we define the continuous function 

n 

M{t):=J2m)\\ (6-20) 

i=l 

then for each t G [0, T*] we may bound 

||S(-,t)||^.(^)<av/M(t) (6.21) 
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for a constant depending on 7^, and n but not on t. 

For if: G [0, T^,] we may argue as in fl6.11l) - fl6.12l) of Proposition l6.2l replacing the temporal 
integration interval [0,T*] with [0,t] to find that 

-/ /to:DS- / Vt, :DH(a,)= / / 5o ■ S + / V 5, ■ S(a,). (6.22) 
Jo Jn Jo Jo Jn Jo -^^ 

We have S(aj(r),r) = and DE{ai{r),r) = for i = l,...,n and r G [0,t], so fl6.22p 
becomes 

To : + So • H. (6.23) 







But by fl6.2ip we may estimate 

- f [ To:DE + So-E< f [ ||S(-, r)||^,(^) (|To| + |So|) 
JO Jn Jo Jn 



< 



r / v^(|To| + |So|). (6.24) 
Jo Jn 



Now we set rj & according to 



7^{s) = —— ==- (6.25) 

2^(1 + ^/M{s)) 

with the constant on the right side of (16.211) : this choice implies that 

r]C^VM <-, and (6.26) 
2 

CVm (^C + = (c + 2a{a + (t)(1 + v^M(i))) < C(l + M), (6.27) 

where we have used Cauchy's inequality, and in the last bound C depends on C^,. We use 
this 7] in Lemma 15.91 to bound 

C,\/M |5o| < r]C,VMCo + C,^ (^C + Ms)ds < j + C{l + M{s))uo{s)ds. (6.28) 

On the other hand, (15.281) of Lemma [5.71 and Cauchy's inequality imply that 

C,^/M \To\ < C{1 + M{s))uo{s)ds. (6.29) 
Now we use flOSjl - flg:^ in fET^ - flg:^ to deduce that 

^* 1^ ■ 5, < ^ £ ^ Co + (cil + Mis)) ^ uo{s)^ ds. (6.30) 

Notice that the right side of (I6.30p is finite since z/q G L°°{[0,T^]; Ai{Q)) and Co has finite 
mass. The bound fl6.30p proves (I6.17P in the case that Fj G C°. If instead Fj G L^, then we 
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let Y^" G C°° be a smooth mollification of Yi so that Y^ — > Yi in L^, and we apply (16.30p to 
K-^ to get 



(is. 



(6.31) 



Then since Y^' — )■ Fj in and Si G L"^, we may send A — )■ on the left side of f l6.3ip : since 
G L°°{[0,T^]; J^{Q)), we may use dominated convergence to pass to the limit A — )■ on 
the right of (I6.3ip . Taking these limits then yields f l6.17p . 

□ 

Next we give a result comparing the measure C to through the use of the "product 
estimate." 



Proposition 6.4. It holds that 







i=l 



5^C> / a7r'J2^'(a,. 



I • |2 



(6.32) 



i=l 



Proof. Proposition 14.41 states that 



Jn 



bV-ijY 



< liminf 

£^0 



Jn 



b IVb^Vs ■ Y 
|loge| 



1/2 



Jn 



|loge| 



1/2 



(6.33) 



for any Y G C^{Q x [0,T,];M2) and G C^{Q x [0,T,]). Note that from f l516|l we may 
compute that 

b l^B^Ve -Yf = Te -.Y ®Y + ge \Yf - |curl B^f \Yf , (6.34) 

from which we deduce, using the structure of V given in Lemma [5.11 and the limits given by 
Lemma [5.51 that 



'^ndib{ai)dj- ■ Y{ 



< 



Jn 



u\Yr + T -.Y ®Y 



1/2 



Jn O! 



1/2 



. (6.35) 



We may apply this result with {Yj,ipj}j^fq a sequence supported in a 7*/j neighborhood 
(with 7* > given by Lemma 15. ip of each of the paths so that 



Yj{ai) — )■ irdi^^^^—^dj' and ipj{ai) — )■ 1 as j — oo. 



(6.36) 



Then since Proposition 16.21 says that Tj = 0, we know that 



Jn 



f n 



Oil 



z/|Fjf + T: Fj-^Fj ^TT^ / fe^(ai)^ as j ^ cx). (6.37) 



i=l 



Also, as i —7- oo. 



ndib{ai)dj' ■ Yj{ai)ds — )■ tt^ / 6 



j=i 



j=i 



(6.38) 
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and 



Jn 



a a 



(6.39) 



The result follows by passing to the limit j — t- oo with these Yj, tpj in fl6.35p . 

□ 

With this and the previous lemma, we can identify u and find that the energy does not 
actually increase by 0(|loge|). 

Proposition 6.5. We have 

n 

F,{v„B,){t) <7iJ2Kaiit))\^oge\ + o{\loge\) for allte [0,^]. (6.40) 

i=l 

Moreover, Co = and for t G [0, T^] we have uoit) = and 

Uiit) = nhi^aiit)) for each i = 1, . . . ,n. (6.41) 
Proof. Integrating the energy evolution equation (14. 7p in time between and t <T^: yields 



|log£:| |loge| 



+ 



Jn 



ah \dtVe 


f + a\dtB,\^ 




\oge\ 





[ [ V{v„B,)-bZ. (6.42) 
Jo Jn 



Passing to the limit e — )■ and using Proposition 16.41 and the well-preparedness assumption 
ffL22|) shows that 

ft ™ I - |2 rt r rt ^ 



ai) ■ 



/ z/(t)- V7r6(a,(0))+ / an'Tb'ia,y-^+ / / Co < / V 27rrf,6(a,)^( 
Ju Jo Jo Jn Jo j^]^ 

(6.43) 

Since G if ^([0, T^,]; fi) and b is smooth, we have that 6 o a-j G if -"^([0, T^,]; M) and hence is 
absolutely continuous on [0,T^,] (cf. Theorems 4.2.2/1 and 4.9.1/1 of |EGj ); this implies that 

-5^7r6(a,(0)) = -^7r6(ai(t))+ / 7rV6(ai) ■ a^. (6.44) 

i=l 1=1 i=l 

By taking the dot product of fl5.39p of Lemma [5.81 with dj, employing Proposition 16.11 for 
Pi = —Uihi, and setting Tj = with the help of Proposition 16. 2[ we find that 

2tx dib{ai) Z {ai) ■ = —2TTdib{ai)Z^{ai) ■ hi 

= di ■ [auidi + z/jVlog6(aj) + (37cdib{ai)dj- + Si] (6.45) 

= aui \dif + Uidi ■ V\ogb{ai) + di ■ Si. 

Plugging (I6.44P and (I6.45P into (I6.43P and using the fact that z/(t) = i'o{t) + ^ i'i{t)6ai{t), we 
find after rearranging that 



n rt r rt IT- 

uoit) + V(z/.(t) - 7r6(a,(t))) + / / Co < a / V 
i=i Jo Jn Jo 



7c^b\a,) 
I a,- 1 I z/,- 



rt n rt n 

+ / V(di- Vlog6(ai))(z/i-7r6(ai))+ / Vd, -S,. (6.46) 
Jo Jo 
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Define the function M{t) = ^"^^ \ai{t)\. We now want to estimate tlie tliree temporal 
integrals on the right of fl6.46p in terms of 1 + M{t). For the first term, we use Lemma [5.61 
to see that TTh{ai) < Vi{t), which allows us bound 



— 



Vi 



{Vi - TTb{ai)) < 2(z/j - Tcb{ai)). 



Hence 



a 



El' 



aA V. 



(6.47) 



(6.48) 



For the second term we use Cauchy's inequality for |aj ■ V log6(ai)| < C(|| V log6||^oo +M) < 
C(l + M), which holds for each i = 1, . . . ,n with C independent of time because b is smooth 
and bounded below by (11. 2p . This yields the bound 



/ V(a,- Vlog6(a.))(^^.-7r6(ai)) < / C(l + M) V(z/, - 7r6(a,)). 
^0 Jo 

To estimate the third term we will use Proposition 16.31 with Yi = hi G L^([0,T* 



(6.49) 



to 



deduce the bound 



fY.a,-S,<\ f I Co+ fc{l + M) ! 
Jo ^ Jo Jn Jo Jn 



Z/Q. 



(6.50) 



Now we sum the estimates (I6.48p . (I6.49p . and (I6.50p and replace in fl6.46p to find that 



/ Mt) + Y.^<t) - nb{a,m + [ 
Jn -^-^ ^ Jo Jn 



< / C(l + M) / z^o + Vi 
Jo Jn 



Vi - nb{ai)) 



(6.51) 



We may view this as the differential inequality 



with 



J^{t) + g{t) < I Q{s)T{s)ds 
Jo 

„ n 

7{t) = / z/o(t) + V(z/i(t) - iib{a,{t))) > 0, 

I = ^ / / Co > 0, and Q{t) = C(l + M(t)) > 0. 
2 Jo Jn 



(6.52) 



(6.53) 



(6.54) 



Note that since G H^{[0,T^]] Q), we have that M and Q are in L^{[0,T^]). We may then 
use Gronwall's inequality on f l6.52p to see that 



^ Q{s)J^{s)ds < (^j^ Q{s)T{s)ds^ exp (^j^ Q{s)ds^ = 0, 



(6.55) 
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and hence that 

J^it) + g{t) < 0. (6.56) 
Since J-", ^ > 0, we immediately see that 



(t) = and Y^{yi{t) - 7r6(ai(t))) = for all t e [0,T,], (6.57) 



the latter of which implies that z/j(t) = 'Kh{ai{t)) for i = 1, . . . ,n since z/j(t) > 'Kh{ai{t)). The 
estimate fl6.40p follows from fl6.57p and the definition of v. We may also deduce from (16.561) 
that that Co = 0. □ 

Since VQ{t)dt and Co vanish, we find that other terms vanish as well. 

Corollary 6.6. We have that Tq = Q, Sq = 0, = 0, and a cnil B^^dtEj- = 0. 

Proof. Since VQ{t)dt = and Co = 0, the vanishing of these terms follows immediately from 
the estimates of Lemmas 15.71 and 15. 9[ □ 

We can use the vanishing of Tq, Sq and T, for z = 1, . . . , n to show that Si = for 
i = 1, . . . ,n. This in turn allows us to compute the dynamical law for Oj and complete the 
proof of Theorem II. 1[ 



Proposition 6.7. The vortex trajectories satisfy ai G C"'^([0, T*]; Q) for i = 1, . . . ,n. More- 
over, they obey the dynamical law 



adi + di(3dj- = -2di{V-^ipo{ai) - ^5^(0^)) - Vlog6(ai) 
= -2di(y^iJo{ai) + V^oiai)) - Vlog6(ai) 

= -2di — - V log b{ai 

V Kai) J 



(6.58) 



Proof. According to Lemma [5.81 and Corollary 16.61 we have that 

n n 

divT = S'o + J] S,{t)5a,^t)dt = Si{t)S,,it)dt. (6.59) 
i=i 1=1 

On the other hand, Lemma 15. 6[ Proposition 16.21 and Corollary 16.61 imply that 

divT = div (to + J2Tiit)^amdt^ = 0- (6-60) 
Equating these two, we find that 

n 

J2S^{t)SaMdt = 0, (6.61) 



i=l 



AA 



and hence Si = for a.e. t G [0, T*]. We use S'j = in the equation for Si given in fl5.39p of 
Lemma I5.8[ and then we substitute in the values of Pi and z/j given in Propositions 16.11 and 
l6.5l to deduce that for a.e. t G [0,T^,], 

ahi + diPaj- = -2diZ^{ai) - Vlog6(ai). (6.62) 

We may solve (16.621) for dj to see that (I1.28P holds. According to Remark 12.51 the vector 
field Z is smooth and bounded, while (II. 2p implies that V log(6) is smooth and bounded; 
then since G C°'^/^, the right side of (ll.28p is a continuous function of t in [0, T^]. Plugging 
(I1.28P into the equation 

ai{t) = ai{0) + / di{s)ds, (6.63) 
Jo 

which follows from the absolute continuity of Oj (again, cf. Theorem 4.9.1/1 of |EGj ). we see 
that ai{t) — aj(0) is the integral of a continuous function and is thus classically differentiable. 
Moreover, (I1.28P implies that ||di||^oo < C so that G C^{[0,T^,]]fl). 

The first equation in (16.581) follows from (16.620 by using Z = Vipo ~ -^o? which is the 
definition of Z given in Lemma [2741 The second equality in (16.580 follows from the definition 
of Xq (I1.12P and equation (ll.lSp . The third follows from Lemma [2.31 

□ 

Although we have only worked on the interval [0, T*], the dynamics can be extended until 
a collision occurs. Indeed, Proposition 16.5! says that 

n 

F,(t;„5,)(r,) <7r5^6(a,(T,))|log£| + o(|log£|). (6.64) 

1=1 

We can then run all of the above analysis again, starting from time t = T^. The only obstacle 
to running this iteration forever is the possibility of a vortex collision or a vortex exiting the 
domain. Hence the maximal time of validity of the theorem is the first occurrence of such a 
collision or an exit under the law (ll.24p . Theorem II. II is proved. 

The proof of Theorem 11.31 can be obtained following the same steps (things are actually 
made simpler by the absence of the gauge B) according to the sketch given in Section 11.41 
Details are left to the reader. 
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